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ABSTRACT 


This thesis conducts an experimental investigation of 
the effects on separation of the closing angle of the stern 
section, the depth of submergence or proximity of the free 
surface, and the speed of a submerged body of revolution. 
The results of this investigation are then compared to the 
results obtained from two existing theoretical criteria for 
the prediction of the separation point on such bodies in 
turbulent flow. The point of separation is observed experi- 
mentally by photographing dye-flow patterns over the stern 
of the model at three speeds from .75 to 1.7? knots; at three 
depths of submergence corresponding to surfaced, near sur- 
face and fully submerged conditions for five different body 
stern shapes increasing in streamlining from a blunt hemis- 
phere to a finely shaped almost pointed section. Equivalent 
bedy profiles to the outside of the turbulent boundary layer 
displacement thickness are computed and potential flow pres- 
sure-velocity distributions based upon the equivalent pro- 
files are calculated and used in evaluating the separation 
criteria. 


It is concluded from the results that initiation of 
separation occurs at a nearly constant body angle on all 
bedies tested regardless of the degree of streamlining. Ex- 
perimental results show no apparent effect of depth of sub- 
mergence on the separation point and the effect of speed 
shews no clear cut pattern. The sensitivity of the variables 
used in the separation criterion on the predictions of the 
separation point are evaluated and recommendations are made 
for further studies and experimentation. 


Thesis Supervisor: Philip Mandel 


Title: Associate Professor of Naval Architecture 
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NOTATION 


Symbol Meaning Units 
C NOon-dimensional coefficients 
E Error functions 
K Matrix functions of body 
ki Matrix functions of prolate spheroid 
L Body length ft. or in. 
Pp One half the perimeter of the body profile ftheorein . 
R Gaussian weighting functions 
Rg Reynold'’s number, along body profile 
Rx Reynold's number, along x-axis 
S Surface area, over body ft.2or in.2 
U Tora city along body ft./sec. 
Wz Free Stream velocity ft./sec. 
¥ Velocity of body ft./sec.or knots 
X Distance along x-axis ft.or in. 
Y Distance along y - axis ft.or in. 
ds Element of arc length along body 
2(x,t) Function of prolate spheroid whose ends coincide 


with body ends and intersects body at x = t 


Acceleration due to gravity 


k Matrix functions of body 

ky Longitudinal virtual mass coefficient of given body 

Kk, (x) Longitudinal virtual mass coefficient of prolate 
spheroid at x = t 

n Reciprocal of exponent in boundary layer vedocity 
profile 

k’ Matrix functions of prolate spheroid 
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NOTATION (Cont'd) 


Symbol Meaning Units 
p Dynamic pressure lb/ft2 
q Stagnation pressure lb/ft? 
r Radius of body profile (calculation in. 


for corrected radius) 


r Distance from location of point in. 
source on x-axis to point on body 
profile (calculation of pressure 
distribution) 


it Location of point source on x-axis, also 
equal to position of Gaussian abscissae on 
X-axis, normalized 


u Component of velocity along x-axis ft/sec 
V Component of velocity along y-axis ft/sec 
Ww Axisymmetric potential function 

x Distance along x-axis, normalized 

y Distance along y-axis, normalized 


Angle between outward normal to body 
and x-axis (calculation of pressure 
distribution) 

Constant in separation criterion 


Angle between normal to body and y-axis or 
between tangent to body and x-axis 


Ratio of §*/s 


Boundary layer thickness in. 
x¢ 

Boundary. layer displacement thickness in. 

Boundary layer momentum thickness in. 


Ratio of body length to body diameter 


< ~ © 


Kinematic viscosity of water ft/sec 
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NOTATION (Cont'd) 


Symbol Meaning Units 
5 Normalized position of Gaussian abscissae 
based on body length of 2 and from -1 to + 1 
IP Mass density of water slugs/ft° 
D Velocity potential function 
yp Stokes stream function 
J Ratio of /§ in calculation of corrected 
body profile 
sO Stream function of a unit source on the x-axis 
Subscript Meaning 
d Subscripts drag coefficient 
f Subscripts coefficient of friction 
i Gaussian abscissae index 


j Gaussian abscissae index 

k Velocity is in knots. 

n Refers to direction of body normal in derivation 
of pressure-velocity distribution 


n Index of iteration for solving equation of 
pressure-velocity distribution 


p Subscripts pressure coefficient 

pl Quantity referred to flat plate 

r Subscripts residual resistance coefficient 
t Subscripts total drag coefficient | 

x Quantity referred to x-axis 
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I - INTRODUCTION 


Separation of turbulent flow with its associated causes 
and effects has for the last half century been a subject of 
controversy in the field of naval architecture and ship de- 
sign. It is in this aspect of fluid flow that the naval 
architect has his greatest lack of knowledge. There is even 
considerable debate and controversy as to just what is known 
and what is not known on the subject of separation, its 
causes and effects. Little has been mentioned in the past 
regarding this aspect of fluid flow in references or texts 
on ship design. What has been stated has been conjecture 
based on limited theoretical calculations, sketchy experi- 
mentation, or recalled experience and observations. 

The causes ‘and effects of separation and the para- 
meters relating to them are complicated and closely linked 
and must be fully understood before any definite design cri- 
teria of value can be established for application to ship 
and submarine forms. Separation of flow on such forms in- 
volves the three dimensional boundary layer and turbulence 
theories which have hardly been more than elaborated be- 
cause, of their associated mathematical difficulties. Even 
the mathematically simpler axialsymmetric flow becomes in- 
volved with the problems of turbulent flow which is not fully 


understood or theoretically formulated for the relatively 
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Simple case of two dimensional flow over a flat plate. 

The ultimate goal is to establish a valid, but simply 
applied, theoretical criterion to predict the separation 
point and a foolproof method of testing such a criterion 
with models to aid naval architects in improving ship de- 
sign. To do the latter the scaling effects involved in 
separation must be firmly established and understood as well 
as understanding the basic phenomenon of separation itself. 

The most important and difficult step in realizing the for- 
mulation of such a criterionis the determination of the sig- 
nificant parameters to be used and the quantitive effect of 
eech on the location of the point of separation. Faired surfaces 
with modest closing ageés apperently retard separation but the 
limiting magnitude of the body closing angle or extent of 
fairing required to minimize separation is unknown. 

The effects of hydrostatic pressure and atmospheric 
pressure on separation in general, and in particular the 
effect on the limiting closing angle have been the subject of 
considerable debate. Saunders (i) *gives some information re- 
lating to these topics but this information has apparently 
never been checked or verified. The effects of the free sur- 
face on separation on a submerged body as it approaches and 
penetrates the free surface is another area of considerable 
debate. Body speed and surface roughness are other possible 
parameters. 


* Numbers in brackets refer to works referenced in Appendix H. 
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The object of this thesis is to investigate the ef- 
fects on separation of varying depth of submergence or proximi- 
ty to the free surface, body speed, and the closing angle of 
the stern of the body. The results are necessarily qualitative 
in nature since the primary tool used in this investigation, 
dye-stream photography,is limited in the exactness to which 
the point of separation can be determined by noting the lo- 
cation of the familiar flow pattern of Separation, i.e., 
abrupt departure of the streamlines of flow from the body 
surface. Originally it was hoped to relate the observed 
movement of the separation point, if any, with the independent 
variation of the parameters listed above to the form drag on 
the body. However, the results of the drag measurements on 
the model were unsatisfactory due to faulty or inadequate 
| a eee and this could not be done. Two separation 
criteria , one developed by Stratford (2) and one aeaucuiies 
by Granville in (3) , are theoretically evaluated and the re- 
sults compared to the observed point of separation from the 
dye-stream photographs. This provides some insight as to 
validity of the criteria and the significance of the para- 
meters used therein. 

It is felt that the results of this thesis are a step, 
although incremental, toward determining what the significant 
parameters in separation of flow are and how they are related 


to the phenomenon of separation. 





I ~ PROCEDURE 


A. General. 


The work presented here was performed in three phases. 
In the first phase model velocity and drag were measured 
experimentally and dye-flow photographs were taken of the 
model stern section for a systematic series of tests. In 
the Bveoad phase theoretical calculations were made to de- 
termine the point of separation on the body which could be 
related to that observed in the photographs. This phase in- 
volved the calculation of the turbulent boundary layer dis- 
placement and momentum thickness, the pressure-velocity dis- 
tribution over the resulting corrected body and the use of 
these quantities in two separation criteria for turbulent 
flow. The third phase was the analysis of the results of 
phases one and two to determine and relate, if possible, the 
effects of the independent variation of speed, depth of sub- 
mergence, and stern shape on the observed separation point 
location on the body as well as to compare and relate the ob- 


served point to the calculated point of separation. 


B. Experimental Phase. 
In this phase experiments to visually determine the 
point of separation on the model and to measure the total 


drag on the model were conducted at the Massachusetts Institute 
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of Technology Ship Model Towing Tank under the supervision of 
the Department of Naval Architecture. Drag and velocity data 
and dye-flow photographs were obtained for the model at three 
speeds, at three a of submergence corresponding to sur- 
faced, near surfaced and submerged conditions, and with five 
varying shaped tail sections which modified the closing angle 


of the model's stern. 


General Description of Equipment. 


Figure I is a sketch of the model showing the variation 
in stern section and the waterlines of the three depths used. 
Figure II shows a view of the towing tank and towing carriage, 
while Figure III shows the carriage control station and in- 
strumentation. 

The model, #4893, was constructed at the David Taylor 
Model Basin, Carderock,Marylandamis a blunt nosed circular 
cylinder consisting of three main parts: a blunt solid ma- 
hogany nose section, a long hollow aluminum parallel middle- 
body, which contains the dye system, and five solid mahogany 
stern sections of varying shape and length. To insure opera- 
tion in turbulent flow a 1/32 inch sieeve _ placed 
around the nose of the model 2.5 inches from the forward end. 
The offsets and other characteristics for the model and stern 
sections are given in Appendix A. Figure IV shows the model 


components. 
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Ficure III - Left to right: carriage contro]. console, 
digital voltmeter and speed measuring equipment, Sandborn 
recorder. 








Figure IV - Model and five stern shapes 





Figure V — Internal view of model parallel middlebody showing 
support plate, dye tubes, dye manifold and brass dye ring. 





The dye system emits dye from 1/16 inch diameter holes 
spaced every 1/4 inch apart around a circular ring located at 
the after end of the parallel middlebody. This system is shown 
in Figure V. Dye was fed to the model from a tank mounted on 
the towing carriage through 1/4 inch plastic tubing. A solu- 
tion of water and potassium permanganate was used for dye. 
Each of the five stern sections can be attached to the mo- 
del immediately aft of the dye ring. When the model was in 
motion dye released from the dye ring would stream over the 
stern section permitting observation of the flow. 

The model was rigidly mounted to the towing carriage 
through a bracket which allowed operation at any' depth within 
the limits of the 18 inch towing street. Figure VI shows 
the model mounted on the carriage. The towing carriage rides 
on an overhead rail and is driven by a variable speed elec- 
tric motor driving a soft rubber wheel in contact with a 
similar side rail. 

| Figure VI also shows the drag measuring dynamometer 
located between the 4:1 ogive towing strut, connected to the 
model. internally, and the 24 inch diameter hollow-cylinder 
which welnconnmcted to the mounting bracket. The dynamometer 
measures the axial force between the carriage and the model 
by means of a differential transformer. 


A length of 32 feet was used for each run. The length 


of run was limited to this length due to a portion of the 
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Figure VI - Model mounted on carriage showing towing dynamometer, 
support bracket, dye tubes and towing strut. 





Figure VII - 4x5 view camera and strobe light arrangement 
at tank window, 
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carriage rail beyond this length being badly distorted and 
thereby cauSing undue carriage vibration. The lenteth of run 
started at the beach end of the tank and extended to the end 
of the first of the 10 foot by 4 foot glass windows located 
in the tank's side wall. A 4x5 view camera was Set up at 
this window along with two high intensity strobe lights which 
were tripped by the carriage striking a contact on the over- 
head rail as it reached the proper position in the window. near 
the end of the run. A white, light reflecting, background was 
installed behind the model because the distance to the far 
side of the tank was too great to achieve satisfactory expo- 
sures. Figure VII shows the photographic set-up that was 
used. Each stern section was marked at 1/2 inch intervals 
along its length as an aid in determining the position of 
separation from the photographs. 

The carriage speed was set by the manual speed control 
dial on the carriage console since the automatic speed con- 
trol was inoperative. Speed data for each run was recorded 
On tape with a Sandbarn Type 150 recorder and type 1300 preampli- 
fier, Drag force was read out on a type 405 Hewlett Packard 
D.C. digital voltmeter with input from the tyye:00 Sandborn 
preamplifier. The voltmeter's purpose was to attenuate the 
high frequency components of the carriage vibrations sensed 
by the dynamometer as drag fluctuations and to allow a volt- 


age proportional only to model drag to be read out. A more 
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detailed description of the towing tank and associated equip- 


ment, model dye system and dynamometer is given in Appendix B. 


Run Procedure 

Two runs were conducted for each condition of stern 
section, depth and speed. On each run a photograph was taken 
of the dye-flow over the stern section while speed and drag 
were recorded simultaneously on the Sandborn recorder and 
digital voltmeter. The tank was allowed to settle between 
each run. Polaroid pictures were taken occasionally to check 
exposure and dye-flow definition. Only five of the dye holes 
in the dye ring were used, two each on the top and the bottom 
and one on the axis of the body. In turbulent flow the dye 
present from the other holes tended to obscure the separation 
point. For laminar flow the use of the entire ring was satis- 
factory. Figure VIII shows both turbulent and laminar flow. 
The unused holes were masked off. All screw holes, etc. on 
the model were filled smooth with a waterproof putty to insure 
as little flow disturbance as possible. 

The model speed was easily obtained from the Sanborn 
trace and a speed calibration curve; however, the drag could 
be obtained from the digital voltmeter only by averaging the 
reading over the run. This was done by recording about ten 
values of the voltmeter reading during the last third of the 


run and averaging. This average reading was then compared 
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to the "zero" drag reading taken prior to each run to determine 
a voltage difference reading proportional to drag. A calibra- 
tion curve was determined by applying a known static load to 

the model and observing the voltmeter reading. Using this curve 


> 
and the wltage difference reading drag force was obtained. 


C. Theoretical Phase. 

In this phase a theoretical position of separation on 
the body was established after correcting the body for the ef- 
fects of bounday layer growth. Throughout this phase it was 
assumed that the flow was turbulent from the nose of the body 
and that the 1/7 power law holds for the velocity distribution 


within the turbulent boundary layers. 


Calculation of Turbulent Boundary Layer Characteristics. 


As indicated in the general discussion of the procedure 
the value of the displacement thickness on the body and mo- 
mentum thickness on the body are necessary. Eckert [4] has 
evaluated the effects of transverse curvature on the turbu- 
lent boundary layer on a circular cylinder in compressible 
flow. These effects are the same as reported by Yu : (5) 


and Landweber (6), which are: 


** The drag measurements were unsatisfactory. Appendix E 


has a more detailed discussion relating to these measurements. 
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(1) Boundary layer thickness and displacement thickness 
are less than on the flat plate. 
(2) Momentum thickness and local and mean skin friction 
coefficients are higher than flat plate values. 
(3) The effects of (1) occur when the ratio of b/r 0.1. 
(4) The friction coefficients deviate from flat plate 
values only when ratio 6/,~ 1.0. 
Eckert develops simple equations which relate the value 
* 
of S, 5” and © on the body to values of the same parameters 
on a flat plate. The assumption is made that the velocity 
profile within the boundary layer conforms to the 1/7th power 


law distribution of U =i and that Blasius‘’ Law relating 
Uoo 


§ 


wall shear stress as inversly proportional to the fourth root 
of the boundary layer thickness are uneffected by the curva- 
ture. The assumption is also made that the boundary layer is 
turbulent from the beginning and it starts from the nose of 
the body with zero thickness. 

The relationships developed for incompressible flow are 
given below and the reader is referred to (4) for the deriva- 


tion details which are straightforward. 


Relationship for boundary layer thickness 


ai eal 
& ai Gs $/3y 
where 


$ pi = 0.38 X (la) 
(Rx) 2°? 


Although (1) is implicit in & it can be easily solved by using 


5 pl 2s a first try and iterating the results. 
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Relationships for Displacement thickness. 


x re. eer 
. 1 V1 + - (a- An \ a 04, ) (2) 


1 





where, * 





* 
= l= 1. (2a) = Spl (2b) 
Amp op 7k 1 é pl 1 +7 
| 
= An C74 1 > Se 
2m Rm) Bnet (ze) = \™ 
ee. O 
Relationships for Momentum thickness. 
9-5 (Qqu pi + . Romer) (3) 


ier. = - Ww = ese = 9 pl (3a) 
PL > 41 TL +2 V1 C742) pl 


AY a eee © (7242) = ye (3b) 
“pl (3741) 7 +1) (27r+1 


* 
The value of S calculated from (2) was applied as a 
correction to the body surface after the manner described in 


(7) as is shown in Figure X taken from (7}. 
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BODY AX/S 


tt = 


Figure X - Equivalent Body 
Profile with Displacement Thickness 
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Thus the corrected body will have a radius given by 
corr = 1+ Cos 8 a (4) 
For use in the calculation of the pressure velocity dis- 
tribution discussed in the next section this radius was norma- 
lized as follows: 


Y corr = corr (4a) 


norm body Length/2 
An IBM 7090 computer program, Program I1, was written to 
evaluateequations (1,2,3,4 & 4a) at 16 points on each of the 5 
bodies for 9 speeds. These 16 points correspond to the loca- 


tion of the 16 ordinates required for the Gaussian 16 point 
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quadrature rule use in the determination of the pressure veloci- 
ty profile. The values of body radius and the angle 6 were 
picked off at each of the 16 points fram full scale drawings of 
the 5 bodies prepared by the authors. A tabulation of this 

data is contained in Appendix A. Details of the computer pro- 


gram are in Appendix C, 


Calculation of Pressure Velocity Distribution 


Potential Flow Pressure and Velocity Distribution was 
calculated according to a method developed by Landweber (8). 
Two types of solutions are given for obtaining the steady, ir- 
rotational, axisymmetric flow of an inviscid, incompressible 
fluid about a body of revolution. The first method involves an 
iterative type solution for Fredholm integral equation of the 
first kind and is not used in this work. The second method ob- 
tains the velocity directly as the exact solution of an inte- 
gral equation of the first kind by an arithmetic procedure. 
Using the second method as outlined in (8] the basic problem 
is to solve for the velocity potential @ which is dependent on 


x, y and satisfies Laplace's equation in cylindrical coordinates: 


=. E)*2 ( ¥ 22) - 0 (5) 
Ox \ Ox/ Oy ony 

where the equation of the normalized body profile is in 
the form: y7 = f (4) (Figure XI shows method of normalizing the 
body) 


and u= - OP , we=- Og 
Ox Oy 








| 

) | 

| | | 

| | 
ze 0 +4 


Figure XI ~ Normalized Body Profile 


There also exists a Stokes stream function Y (x,y) 


which satisfied the equation: 


O2Y 4 D2Yu 1 DY (6) 
sf 58° 4 oF 


For uniform flow of velocity (Z parallel to the x-axis 
G= -U[x and VY = = 4Uly? at infinity 
The following boundary condition on the body must also be 


e. = = V cos& 


; aw) S ‘ 
Where V = velocity of the body and <= the angle between the 


satisfied: 


outward normal to the body and the x-axis. 


Derived from Green's Theorgnthe equation: 


Jomas- fe @ as (7) 


1g 





where @ and w are any functions harmonic in the region exterior 
to a given body and vanishing at infinity, is used to obtain 

an integral equation for axisymmetric flow about a body of 
revolution. 


Substitution of: dS = 27 y ds and 


ag = - sin s into (7) gives: 
dn 
p Pp 
[ y®@ dWuas 7 y wsin ¥ ds (8) 
O dn ps 


y = ordinate of meridian section of the body which is 
moving in the negative x — direction with unit velocity. 

ds= element of arc length along body 

i= angle of the tangent to the body with x =~ axis 

2P= perimeter of a meridian section 

With Was an axisymmetric potential function and 


Y (x,y) the corresponding stream function 


y dw= dy 
dn ds 
and Pp P Pp 
y® dw ds ™ gy -[-Y 28 as 
5 dn O O ds 





If (f is total velocity along a body when a stream of 
unit velocity is superposed in the positive x - direction 


([ =-dg + cos7x 
ds 


Since dx = ds cos and dy = ds sin X equation (8) can 


be written as: 
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P /P | P 
‘i Uy ase/ (Y dx - yw dy) - py (9) 
O O Y 


but y= “ = ~ of since w and Y are corresponding potential 
Zs y 


and stream functions and #/dx - ywdy is an exact ‘differen- 


tial defining a function {2 (x,y) such that 


a 2 a 

OF oY 

From this and y Ow = OY we obtain 
OY Ox 

O72 +97Q +1 aN 

eS) wayne. Be: y O y 


which is identical to (6) and is the equation satisfied by 
Stokes stream function. It can also be verified that w and V 
are corresponding axisymmetric potential and stream functions 


which satisfy the equation: o¥ =-y O¢ , d¥n=y $s 


Ox OY dy 
Writing (9) in terms of () 
P | P (10) 
wr Of2 as = (£2 - @ ae 
5 OF . O¥ 


If fLis the stream function of a point source of unit 


strength situated at x = t @the x - axis 





SU = - 1 + x-t » A= \(x-t)? + y2 | (11) 
Se 
Then OS = y" , and if y vanishes at both limits, 
Ox st 
Pp 
(Q- g oN ) = 2 
% O 
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(11) now becomes: 


{# (x)y2 (x) ds = 1 (12) 
at 
Oo 


which is an integral equation of the first kind in which the 
unknown function is (4 (x) and the kernel is Va oe When CZ 
is determined the pressure distribution is given from 


Bernoulli’s equation. 


P= 1 - (v2 w?2) = 4-0" 


where q is the stagnation pressure, 
Landweber uses an iteration procedure to solve (12) 
as follows: 
If the iterations are expressed in terms of error functions 


n (t) at the point x = t , first obtain 


1 (t) from yn (t) = 1 -[- (Ln (x) oa (x) ds 
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by using the first approximation //, (t) = (1+k) cos ¥ (x) 
where k, = longitudinal virtual mass coefficient. Then suc- 


cessive iteration Bo, Eg «2106. E, are obtained from 


- 
Eo) (t) = EL (t) - af En_(x)_y 2(x) dx 
oe 


where “o » Ky are nose and tail absciss#of the body.and 
finally Uf n4;(t) is obtained from 


n 
Uni, (t) = UY, (tt) + cos ¥ (t) 2 By Ct) (13) 


p hea 
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To evalyate integrals of the form 


x 
f =) = 
Xo St 


use k (x,t) = # (x) (14) 


| (x—t)2 4+ F (x) B72 
Kk’ (%t) = (%, t) (15) 
|< ay. + @(X) B/2 


where y? = Ff () is the equation of the given body profile 
and y4 = g (x, t)) is the equation of the prolate spheroid 
whose endS cOincide with the ends of the given body, and 


which intersects the body at x = t, i.e., 


g(x, t) =f (t) _(x-x0) (x, -x) (16) 
t= Xo xy — t 


| | 
For the spheroid the length - diameter ratio A is given 


by a2 = (t - xo) (x - t) (17) 
f~ t 


and from this the virtual mass coefficient k, (t) can be ob- 


tained from: 


k, = Aln(A+/22 -1) - rA~ 1 
M2774 -1 - Pin(A+ 4 1) (18) 


E, (t)is then obtained from 


i x; 
E, (t) = 1 — 1+k,/ [k (x,t) -kK’ (x,t)| dx - 14+ k& (19) 
i = [2 ) ) ] TK GY 


and finally ie 
E,4, (td) © E At) E,(t) -$ | K(X, t) E, (x)-Bp (t) |ax 
Yo (20) 
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In using Landweber'’s method the velocity is obtained 
at 16 points along a body, from t = -1 to t = + 1, correspond- 
ing to the Gaussian values CAL for the 16-point guadrature 
rule. In applying the rule the integrands in (19) and (20) 
are evaluated at the 16 Gaussian abscissae Xj = Ke for each 
of the 16 values of tj=%$4 . Gaussian weighting functions 
are used to obtain K j i 
(X;, ti) which are used in evaluating (19) and (20) .U,(t) is 


/ 
= Rj k (Xj, ti) and Kj 4 = Rjk 


then given by (13) and Pla by Bernoulli's equation in the form 
we -1-U, - 

In the application of Landweber'’s method of determining 
the pressure-velocity distribution on a body of revolution to 
the models used in this thesis the corrected body profile as 
calculated in the previous section was used. The function 
Y? = (x) for the normalized corrected body profile was not 
determined because even if it is determined the end result is’ 
still a point-wise distribution of U/U,. Instead of determin- 
ing the polynomial describing the body profile the values of 
the normalized corrected body radius at each Gaussian ordinant 
computed in Program II were used. The angle, 0, was again de- 
termined from full scale drawings of the model and stern sections. 

An IBM 7090 computer program, Program III, was written to 
determine the préésure-velocity distribution on the body cor- 
rected for boundary layer effects using Landweber's method. The 


program calculates C U/Uloo, absolute velocity and gage pres- 


p’ 
sure on the body at 16 points on each of 5 bodies for 3 different 
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speeds at each of 3 different depths of submergence. The 9 
speeds used are the same as those uSed in Program II and are 
the average of the speeds listed for each run and its dupli- 
cate, i.e., Run 1 and Run 1A, in the tabulation of experimen- 
tal data in Appendix D. Thus the outputs of the two programs 
are a compatible package for application in the separation 
criterion to be evaluated. This output is in turn compati- 
ble with the dye~-flow photographs since the same conditions 
apply. A test solution for Program III was run on the ex- 


ample worked out in [7]. 


Separation Criteria . 

Stratford (2) has developed a rapid method for the pre- 
diction of flow separation in turbulent flow from an approxi- 
mate solution of the equations of motion. The equations are 
integrated by a modified “inner and outer layer solution" 
technique developed by Stratford in 1954 for laminar flow 
boundary layers. The final solution is given by equations 
(21) through (24) which relate X, Cp, and dCp/dX at the sepa- 
ration position and may be applied to any given pressure dis- 
tribution. The algebra and procedure for obtaining the sepa- 
ration condition are clearly detailed in {2} and only the re- 
Sults will be given here. 


The resulting implicit equation for the separation po- 


20 








Sition, X, is i - 1 
£ 160 (21) 
(2cp)4 2) (x ach = 1.06 G (107° Rs) 
w n-2 
for Cp = 
n+l 


Where: Ro is the Reynolds number based on value of distance 


X and freestream velocity, Uoo. 


; n= log 10 R, (22) 
and &= 0.66 for d“p < 0 (23) 
g@xe 
f = 0.73 for d“p > 0 (24) 
dx 


with the resultant value of ap being that immediately prior 
dx 


to separation. This formulation has a range of uncertainty 
of 10% due to its inability to fully account for d“p/dXx in 
the case of turbulent boundary layers. MWMowever, the rapidity 
of the method makes this range of uncertainly acceptable in 
this case. Granville in (3) has brought forward another even 
Simpler separation criterion It was apparently originally 
developed by Bui in 1931 at Zurich. The criterion states 


that separation will occur if the following relation holds: 


®,UVoo , du, L < -0,006 (25) 
u L dX Uoo ~ 


Where: 9 = momentum thickness on the body at separation 


point (Calculated using equation (3) ) and L = body length. 
) 
Since L and Uoo are constants for each body at a parti- 


cular speed 


>) 
an 4 lee Uoo 








Equation (25) can thus be written 
©; , Yoo . (gas -)= - 0,006 (26) 


; ary 


Unfortunately the information regarding this criterionis 





limited to that obtained by the writers from Course 13.012 in 
the Department of Naval Architecture at M.I.T. and the meager 
yet reliable reference to it in (3} . A range of uncer- 
tainty for the criterionis not known. It is intended to be 
used as a criterionto be used in scaling separation from 
model to ship seale. 

| With the velocity profiles, pressure coefficients, and 
momentum thickness values on the body determined as outlined 
in the previous sections it is possible to evaluate both of 
these criteria and see how they will compare with the experi- 
mentally observed points of separation. This was done by 
hand calculation for a few cases. Sample calculations are given 
in Appendix F. The results of such calcuations are in the 


RESULTS section of this thesis. 
D. Analysis Phase. 


The data obtained from the experimental phase relating to 
drag was analyzed by utilizing an IBM 7090 computer progran, 
Program I, to reduce the data to coefficient form. The details 
of this computer program are given in Appendix C. The analysis 
of the experimental drag data contained in Table D-I of Appen- 
dix D showed the data to be unreliable and faulty. Appendix E 


contains some of the results obtained from the computer program 
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and a discussion relating to these reSuits. The second step 
in analyzing the experimental data was to study the dye-flow 
photographs and determine from them the point of separation 
and the body angles with respect to the body axis at the point 
of separation. The point where the dye-flow appeared to a- 
bruptly depart from the body profile was taken to be that of 
the point of separation. Only the flow along the keel of the 
model was used in determining the location of the point of 
separation because it was felt that the presence of the towing 
strut forward would influence the flow over the top of the 
model. The authors each made independent determinations of 
the separation position on each run. Differences which arose 
in interpreting the location of the separation point were 

then compromised. The determinations were also made prior to 
the calculation of the separation position from the two sepa- 
ration criteria so that the location of these theoretical 
points would not influence the interpretation of the photographs. 
In this manner it was felt that as high as possible degree of 
precision would be obtained in interpreting the photographs. 
There was, as was to be expected, considerable judgment in- 
volved in establishing the separation point in several of the 
Pictures. A sampling of typical photographs from which the 
separation point was determined is shown in Figure XII. The 
location of the observed separation ‘point for each run is con- 


tained in Table D2 in Appendix D. The angle between the body 
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Run 7 Run 8 Run 9A 


A. Sequence of Runs on Body 3. 








Run 35A 
Body 4 


Run 20A 
Body 2 


R, Runs on Cther Bodies. 


Figure XII - Typical Dye-flow Photographs 
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surface and the body axis at the point of separation was de- 
termined graphically by plotting the separation point loca- 
tion on full scale drawings of the bodies. The movements of 
the point of separation were analyzed to see what, if any, 
correlation could be made between these movements and changes 
of speed, depth of submergence or body shape. The angle 

at which separation occurred .was similarly analyzed. 

The third step in the analysis was to calculate the 
point of separation by each of the two separation criteria 
using the ne of computer programs I and II. Samples 
of the output from these programs are contained in Appen- 
dix G. Sample calculations for each of the criteria are 
given in Appendix F. During these calculations the criti- 
cality of the various parameters in the criteria were 
noted. The results of these calculations were compared 
to the observed point of separation noting location and 


angle of body surface at the point of predicted separation. 
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III RESULTS 


The results of the procedures described in the pre- 
vious sections are presented in this section. These results 
are discussed and conclusions are drawn from them in the 
next section of the thesis. 

The results are presented primarily in graphical 
form. Figures XIII through XVII show the relationship be- 
tween the pressure-velocity distribution over the body and the 
location of the observed separation points. Figures XVIII 
through XXII show the location of the experimental and calcu- 
lated points of separation. The body profiles corrected for 
boundary layer effects are also shown here. Table I con- 
tains a breakdown of the experimentally observed separation 
points on each body by depth of submergence, and speed. The 


body angle at separation is also shown. 
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TABLE I 


Results of Studies of Dye Flow Photographs and Separation Criteria 
BODY #1 X / 
L 


Code Angle Points (Speed-Depth) 
— 348.90 @ LM, 24 LD, 2@ MD, 2 @LS 


B - 32°  ,970 3 @ MM, 2 @ HM, 1 @ MS 
f= 429 #4978 2@HD, 2@ HS, 2 @ MS 
E - 60 945 BURI CRITERION 
F - 47° 982 STRATFORD CRITERION 
BODY #2 y 
Code Angle Jy, Points (Speed-Depth) 
Bo 967 1@LD, 1@LXx °° 
Ee —- 30° «4.976 2@MM, 1@HM, 1@ LD, 2 @LS 
Cc - 36° .983 2@MD, 1 @HD, 1 @ MS, 1 @ LM 
D - 42° .989 2@uHS, 1 @ MS 
E - 13° .930 BURI CRITERION 
F - 389 . 986 STRATFORD CRITERION 
BODY #3 x 
Code Angle ~/L Points (Speed - Depth) 
| =e ~ ,967 2@MS, 2@LS, I @uHs 
B - 285 975 2@LM, 2 @ MM, 1 @ HM, 1 @ HS 
= —- 33 982 2@HD, 2 @MD; 2@LD 1 @ HM 
E - 10°  .905 BURI CRITERION 
F - 46° ,992 STRATFORD CRITERION 
BODY #4 x 
Code Angle /Y, Points (Speed - Depth) 
- 2 983 2@00), 1@wDb °°  — 
B —- 38° 990 2@LS, 2 @MS, 1 @ HD, 1 @ MD, 1 @ LM 
C - 49° .,996 2@HM, 2 @ MM, 2 @ HS, 1 @ LM 
= - 6° .881 BURI CRITERION 
F - 36° . 989 STRATFORD CRITERION 
BODY #5 x 
Code Angle fu Points (Speed=-Depth) 
A - 330 9978 All speeds at deep depth 
B - 30° # .988 All speeds at medium depth 
C - 60° ,998 All speeds at surface 
me —- 5° 831 BURI CRITERION 
F - 23° #£«.979 STRATFORD CRITERION 


Separation points are plotted on body profiles in 
figures XVIII through XXII. 


Point designations are the same as outlined in Table D-II 
of Appendix D; e.g. 2 @ HM means two points were plotted at 
that location and angle on the body corresponding to two runs 
at the high speed —- medium depth condition. 
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IV DISCUSSION OF RESULTS 


Accuracy of Observations 


The accuracy achieved in determining the location of 
the separation points from the dye-flow photographs was 
limited by the ability of the authors to pick off the point 
where the flow abruptly departed from the body surface and 
in their ability to interpolate between the 4 inch markings 
on the bodies. In some cases the point of separation was 
partially obscured by excessive amounts of dye while in 
Other cases, especially at the higher speeds, a portion of 
the dye appeared to be carried beyond the separation point 
by flow outside the boundary layer. In the surface condition 
there was some evidence of dye being drawn up into the area 
aft of the body toward the free surface, The accuracy in 
observing the point of separation in turbulent flow can not 
be equal to that obtainable for laminar flow, but improvement 
in this accuracy could be gained by using a larger model and/or 
using a more sophisticated means of photographing the flow 
such as telephoto lens or underwater motion pictures to ob- 


tain greater definition of the flow about the body. 


Effect of Body Closing Angle 


In examining the positions of the observed points of 
separation on the various bodies as shown in Figures XVIII 


through XXII it was noted that the points occurred in groups 
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located at three or four positions on the bodies. Table I 
shows the details for these groupings. The center of each 
group was used as the measuring point of the body angles. The 
degree of accuracy obtained in measuring these angles is a 
function of the accuracy achieved in the initial location of 
the point of separation already discussed and the errors en-~ 
volved in measuring the angle from the full scale drawings. 
The two factors result ina possible error of = 3 degrees. 

The initial group of separation points on all five 
bodies tested occurs at a body angle of between 23 and 28 
degrees, averaging about 25 degrees. There is no apparent 
orderly correlation among the five bodies to account for this 
variation except that it is due to experimental error. This 
results in the conclusion that the initial point of separa- 
tion occurs at approximately the same angle for all bodies 
regardless of the amount of streamlining. A critical angle 
of 25 degrees thus appears to be established for initiation of 
separation. The rate at which the body angle is changing does 
not appear to effect this critical angle because bodies 1 and 95 
which have vastly different rates of change of body angle have 
nearly the same critical angle. A systematic check over a 
greater range of body shapes, speeds and depths of submergence 
is warranted to check this conclusion. 

The maximum value of the angles in the region over 


which the separation points occurred varied from about 42 degrees 
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to 60 degrees showing that the observed point of separation 
did fluctuate. At these higher angles, however, the speeds 
were generally toward the high end of the speed range and the 
fluctuations may be caused by dye being carried aft by flow 
outside the boundary layer, especially on body 5S where the 
point of separation was most difficult to determine. Be- 
cause Of this possible effect these angles are not con- 
sidered as significant as the initial separation angles. 
The above Soatent:. notwithstanding, these higher angles do 
point up that under certain circumstances separation of 
flow may not occur until the critical angle has been ex- 


ceeded by a considerable amount. 


Evaluation of Separation Criteria 

Buri's theoretical separation criterion equation (26), 
consistently predicts the separation point well forward of 
both the experimentally observed points and the point of 
separation established by Stratford's criterion, equation 
(21). For Buri's criterionthe body angle values all fall 
between 5 and 13 degrees el no apparent consistency of 
vayiation, except that when body 1 is eliminated the angles 
decrease from body 2 to body 5, thus indicating that the 
Peadr ation point moves forward with increase in streamlining. 
This is inconsistent with the sais Bp eel results and the 


general concepts relating to separation which maintain that 


streamlining should decrease the extent of the separation zone, 
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not increase it. This criterionis very sensitive to slope of 
the curve of U/Uoe versus X/t, as plotted in Figure F-I, Appen- 
dix F, and relatively insensitive to model speed that enters 
through 9 in (26). The criterion predicts separation to occur 
gust as the value av/vod) /a (X/;)) becomes negative. A slight 
error in plotting this curve or in using it to solve (26) 
could result in a larger error in the separation point; how- 
ever, the magnitude of the possible error, 20.2 inch, is not 
sufficient to improve the agreement between Buri and the ob- 
served points. 

On the other hand, Stratford's criteriongives much 
closer agreement with the experimentally observed points of 
separation, as it consistently predicts separation either with- 
in or slightly aft of the observed region of separation. The 
values of the body angles at the predicted separation points 
vary from 23 to 47 degrees. Again there is no apparent con- 
Sistency of this variation. Depending on the shape of the curve 
of Cp versus X, as plotted in Figure F-II, Appendix F, this 
criterionis sensitive both to the value of Cp and the value 
of d(C vax and also insensitive to model speed that enters 
through R_ and n in (21). Where the slope is steep the value 
of Cp is the controlling factor while for a shallower slope 
both Cp and d(Cp/Ydx control. The criterionessentially pre- 
dicts separation to occur very shortly after Cp becomes posi- 


tive. Again errors in plotting or interpreting Figure F-II 
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woulda not effect the predicted point of separation by more than 
=.0.1 inch. 

It can be concluded that Buri's criterion did not give a 
b@eistactory comparison to the observed point of separation 
while Stratford's agreed much better although it generally 
indicated that the body angle at separation was of a larger 
yvatue than that observed. It is to be kept in mind, however, 
that these criteria were analyzed on the basis of the calcu- 
lated profiles in Figures AIII through AVII which were de- 
rived tov modified potential flow. it is therefore recom- 
mended that an experimental determination of velocity-pres- 


sure protiles be made on the bodies and that these be applied 


to the separation criteria before a final judgment is made. 


Effect of Vroxinity of ree Surface 





Table I yields no apparent consistent variation of the 
experimental points of Separation on the various bodies with 
a change in oporating depth. 

Computer Program IIT calculates gage pressure or dy- 
namic pressure plus hydrostatic pressure due to depth of sub- 


my 


mergence on each body. The variation ih pressure at the keel 
oS from 2bout 19 lbs. per square foot at thé surfaced condi.— 

eon to about “95 lbs. por Square .feot at the deep depth whach 
Shoutd be scufticient to show the effect of change in depth on 


the observed tocation of the separation point if in fact such 


an effect on separation does exist. Loth of tne separation 
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criteria are completely independent of depth since it does not 
enter into any of the variables contained in equations (21) or 
(26). It is concluded from the experimental results that hydro- 
static pressure over the range investigated has no recognized 
consistent effect on the location of the separation point. 

It is recommended that further experimentation be 
conducted with larger variations in depth. Perhaps the use 
of a water tunnel or propeller tunnel would be better than 
a towing tank for this type of experiments since a large 
static pressure could be applied to the model in the test 
section. More typical ship shapes rather than body of re- 
volution forms should be used in conducting dye~-flow experi- 
ments in the air-water interface to show the effects of 


draft on separation over a given form. 


Effect of Speed 


Bodies I through 4 show values of the body angle at 


the separation point between 23 and 28 degrees for the low 
speed range, between 30 and 37 degrees for the medium speeds 
and 33 to 60 degrees for the high speeds. Because of the dif- 
ficulty in locating the separation point for body 5, no rela- 
tion between speed and separation point could be seen for that 
body. The results for bodies 1 through 4, however, show that 
the separation point appears to move aft on the body as speed 


increases. The reason for this effect could be due to the 
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fact that dye was being carried back along the body outside 
the boundary layer and misinterpretation of the photographs 
resulted. At the speeds used the Reynolds' Numbers were from 
3.5 to 7x 10°. This is within the doubtful transition region 
between laminar and turbulent flow and the possibility exists 
that the lower speeds were at the laminar end of the transi- 
tion, although an examination of the photos in Figure XII 
indicates a typical turbulent flow pattern. As speed in- 
creased the mere fully turbulent flow region was approached 
with the added increased effect of the cmap wire with speed 
thus causing the point of separation to move aft on the body, 
as it does ina typical transition from laminar to turbulent 
flow. The two theoretical separation criteria are insensi- 
tive to changes in speed in the range investigated as has 
been previously discussed. 

The recommendation here is obvious; the experiments 
must be conducted at higher speeds to further insure turbu- 
lent flow. In the present experiments the speeds obtainable 
were directly related to the ability to observe the dye flow, 
the dye being dissipated with increased speed and almost 
wholly indistinguishable beyond 1.75 knots. The speed limi- 
tation was due in part also to large carriage vibrations at 
higher speeds than those used. With a more refined method of 
turbulence stimulation and dye-flow obsbkvattion techniques 


higher speeds could be attained leading to perhaps more con- 
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clusive and plausible results relating the effect of changing 


speed on the separation point than those herein achieved. 


Corrected Body Profiles and Pressure-Velocity Profiles 


Equivalent or corrected body profiles shown in Figures 
XVII through XXII are the results of Computer Program Il. 
From the sample of these results given in Table G-I in 
Appendix G it can be noted that there is only a very slight 
variation in corrected profile on a given body with varia-~ 
tion velocity. The corrected body profile radii decreased 
slightly with increase in speed. The variation is considered 
to be beyond the accuracy of the model construction and the 
full scale drawings so the corrected body profile was plotted 
the same for all velocities on a given body. In Figures XVIII 
through XXII the corrected body is plotted beyond the separa- 
tion point however, this in fact does not exist since this 
profile is based upon modified potential flow theory which 
does not hold after separation occurs. 

The pressure-velocity profiles for the bodies shown 
in Figures XIII through XVII are the results of Computer Pro- 
gram III. From the sample of these results given in Table G-II 
in Appendix G only a slight variation in Cp and Nes occurs 
with variation of Uoo or model velocity. This is a result of 
the aforementioned slight variation in body profile with model 
velocity since the output of Program II, Y — NORM, is the in- 


put to Program III. Again this slight variation is considered 
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beyond the accuracy of the model drawings and certainly is 
not detectable when plotted at the scales used in Figure XIII 
through XVII. Thus the profiles appear the same for all ve- 
locities on a given body. 

Since neither depth of submergence nor hydrostatic 
pressure enter into the calculations of the corrected body 
profile or pressure-velocity profiles none of the figures 
show a variation with depth. 

In summary, the experimental results indicate that 
separation appears to begin at a critical angle of about 
29 degrees irrespective of body streamlining. The correla- 
tion of separation with speed and depth is not clearly seen 
since the speed and depth range should be larger for more 
conclusive results. The results as pointed out in the in- 
troduction are qualitative in nature and are not be be con- 
strued as quantitative evidence. Such evidence must be pro- 
vided in later works. The primary recommendation is that 
further agitied be carried out at greater speeds with im- 
provements in the technique of observing the separation 
point. Stratford's separation criterion, which is more ri- 
gorous than Buri'’s appears to give good agreement with the 
experimental observations, but it should be further evaluated 
by using experimentally determined pressure-velocity dis- 


tributions. 
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V — CONCLUSIONS 


The conclusions which were arrived at in IV are summa- 
rized as follows: 

1. The observed position of separation was determined 
within an accuracy of = 1/2 inch on the body and the body 
angle at this position was eneebe determined to 2 3° 

2. Initial separation occurs at body angles of be- 
tween 23° and 28° regardless of the amount of body streamlining. 

3. Buri'ts criterion for separation consistently predicts 
separation well forward of the observed separation point and 
thus does not compare favorably with the observed results. 

4. Stratford's criterionevaluated from the calcu- 
lated pressure-velocity distributions over the bodies gives 
reasonably close agreement with the location of the observed 
separation point. 

5. Results of analyzing the separation position lo- 
cation show that it moves aft with increase in speed but this 
may be due to the fact that at the speeds investigated the 
Reynolds numbers were in the transition region between laminar 
and turbulent flow. 

6. The theoretical criteria evaluated are both in- 
sensitive to body speed changes and are independent of the 


depth of submergence of the body. 
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7. The observed separation point does not move in a 
consistent manner on the body with a variation in depth of 
submergence of the body. | 

8. The corrections applied to the body surface to 
account for boundary layer effects showed little change with 
change in body speed for a given body. 

9. The pressure-velocity profiles over the bodies 
calculated from modified potential flow showd little change 


with change in body speed for a given body. 
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VI — RECOMMENDATIONS 


The recommendations made in IV are summarized as 
follows: 

1. Further experimentation should employ a larger 
model if possible and an improved means of observing the 
flow about the separation point. 

2. <A greater range of speeds and depths should be 
investigated in future work to increase the reliability of 
the data obtained. Higher speeds consistent with the a- 
bility to observe the dye-flow are necessary to insure that 
flow is definitely in the turbulent region. 

3. If higher speeds can not be attained improve=- 
ment should be made in the method of inducing turbulence 
through stimulation. 

4. -If dye studies are continued a minimum of dye 
should be used to prevent dye outside the boundary layer 
from obscuring the observations. 

5. Pressure measurements should be made in and 
ahead of the separation zone to experimentally determine 
the pressure-velocity distribution over the bodies and also 
serve as another method of locating the separation point. 
This could be compared with the results of this work and 


used to re-evaluate the separation criteria on the basis of 
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the actual pressure-velocity profiles rather than on calcu- 
lated profiles. This appears to be the next logical step 
as the initial determination of the separation zone has been 
made in this work. 

6. Recommendations relating to drag force measure- 


ments are included in Appendix E. 
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APPENDIX A 
OFFSETS AND OTHER CHARACTERISTICS FOR MODEL AND STERN SECTIONS 


TABLE A-I -— OFFSETS 





o7 


Body: Nose Section Stern 1 Stern 2 Stern 3 Stern 4 Stern 5 
Station Spacing: 0.50" 0.50" 0.50" 0.50" 0.50" 0.50" 
Station 0O at: Nose Tail Tail Tail Tail Tail 
Station: Radius(Inches) Radius Radius Radius Radius Radius 
0 0.00 0.00 0.00 0.00 0.00 
0.5 0.43 0.80 0.65 0.63 0.32 
1 0.68 a 1.05 0.87 0.85 0.47 
£.5 0.83 x 1323 1.03 0.98 0.57 
2 0.94 S Nes yf 1.16 1.09 0.69 
3 Lea Z ~ 1.58 1.36 L viet 0.87 
4 1.27 a) 1.72 Lwo2 1.41 1.00 
5 1.37 1.84 1.66 1.53 1,12 
6 1.47 1.93 1.77 1.63 1.23 
7 1.54 1.98 1.84 1.72 L<3s 
7.26 ae 3 2.00 = = == 
8 1.62 a 1.91 1.78 1.42 
9 1.67 < 1.95 1.84 1.50 
10 1.72 = 1.98 1.89 1.57 
1l 1.76 Tt 1.99 1.93 1.63 
11.30 aa 0 2.00 -- -= 
12 1.80 1.96 1.69 
13 1.83 ve 1.98 1.74 
14 1.87 Ly 1.99 1.79 
14.78 — a 2.00 -= 
15 1.90 ln 1.83 
16 1.93 = 1.87 
17 1.95 Li 1.90 
18 1.96 < 1.92 
19 1.97 1.94 
..20 1.99 1.96 
21 2.00 1.98 
22 ie 1.99 
22.78 os 2.00 





TABLE A-Ii- OTHER CHARACTERISTICS 


Nose/Tail Length Wetted Surface Wetted Surface 


Radius (inches) Depths 1 and 2 Depth 3 

(inches) (sq. ft.) (sq. ft.) 
Nose 0.43 10.5 0.726 0.640 
Cylindrical Section ae 19.0 1.653 1.31 
Stern l 2.00 2.00 0.175 0.158 
Stern 2 1.47 3.63 0.245 0.219 
Stern 3 0.98 9.65 0.387 0.340 
Stern 4 0.99 7.39 0.513 0.454 
Stern 5 0.33 11.39 0.741 0.664 
Body l =-- 31.50 2.904 2.109 
Body 2 a= 33.13 2.624 2.170 
Body 3 ~= 30.15 2.766 2.291 
Body 4 a 36.89 2.892 2.405 
Body 5 = 40.89 3.120 2.615 


Body means total of nose, cylindrical section, and stern shape. 
Wetted surfaces submerged (depths 1 and 2) were calculated 

using ten-ordinate Simpson's rule with half-stations on body perimeters. 
Spacings of 1 inch were used on nose and stern 5, 3/4 inch on 
stern 4, 1/2 inch on sterns 2 and 3. Stern 1 was calculated 
directly since it was a hemisphere. 
Surfaced wetted surfaces (Depth 3) were calculated by computing 

non-wetted surface using trapezoidal rule on non-wetted arc lengths 

and subtracting from submerged wetted surface values. 


These calculations are not shown in this report. 


58 





TABLE A-III 


GAUSSIAN STATION,RADIUS, AND BODY ANGLE 
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Gaussian X Bod | Body 2 Body 3 
Station L xX Radius Angle X Radius Angle X Radius ngle 
-.990 .005 0.17 0.37 49.5 0.18 0.37 49.4 0.19 0.40 47.3 
-.944 .028 0.87 0.88 26.0 0.92 0.91 25.0 0.97 0.94 23.9 
=-.866 .067 2.12 1.30 14.2 2.23 1.31 12.7 2.36 1.37 11.7 
-,756 .122 3.85 1.60 7.5 4.05 1.63 7.1 4.30 1.67 6.3 
-.618 .191 6.02 1.81 4.5 6.33 1.84 4.4 6.72 1.87 3.6 
-.458 .271 8.54 1.96 2.6 8.98 1.98 1.8 9.52 1.97 1.9 
-,282 .359 11.31 2.00 0.0 11.90 2.00 0.0 12,62 2.00 0.0 
-,.096 .452 14.25 2.00 0.0 14.99 2.00 0.0 15.90 2.00 0.0 
+,.096 .548 17.25 2.00 0.0 18.14 2.00 0.0 19.24 2.00 0.0 
+.282 .641 20.19 2.00 0.0 21.23 2.00 0.0 22.52 2.00 0.0 
+#.458 .729 22.96 2.00 0.0 24.15 2.00 0.0 25.62 2.00 0.0 
+.618 .809 25.48 2.00 0.0 26.80 2.00 0.0 28.43 2.00 0.0 
+.756 .878 27.65 2.00 0.0 29.08 2.00 0.0 30.85 1.94 5.1 
+.866 .933 29.38 2.00 0.0 30.90 1.78 13.1 32.79 1.62 15.5 
+.944 .972 30.63 1.63 35.0 32.21 1.33 27.8 34.18 1.15 24.5 
+.990 .995 31.33 0.78 67.2 32.95 0.69 59.6 34.96 0.97 51,5 
Gaussian X Body 4 Body 5 

Station L X Radius Angle X Radius Angle 

-.990 .005 0.20 0.41 47.1 0.21 0.42 46.5 

-.944 .028 1.02 0.97 22.7 1.13 1.00 21.6 Wi 

-.866 .067 2.48 1.38 11.5 2.75 1.42 10.0 7. 

-.756 .122 4.51 1.67 5.%7 5.00 1.72 ae X is inches from 
-.618 ,191 7.05 1.88 3.4 7.81 1.91 2.0 nose of body. 
-.458 .271 10.00 1.99 1.0 11.08 2.00 0.0 Radius is in 
-.282 .359 13,25 2.00 0.0 14.69 2.00 0.0 inches. 

-.096 .452 16.69 2.00 0.0 18.50 2.00 0.0 Angle is in 
+.096 .548 20.20 2.00 0.0 22.39 2.00 0.0 degrees. 
+.282 '.641 23.64 2.00 0.0 26.20 2.00 0.0 

+.458 .729 26.89 2.00 0.0 29.81 1.98 1.0 

+.618 .809 29.84 1.99 1.0 33.07 1.85 4.1 

+.756 ..878 32.38 1.84 6.2 35.89 1.597 8.3 

*+.866 .933 34.41 1.52 12.0 38.14 1.18 -’12.9 

+.944 .972 35.87 1.10 22.0 39.76 0.73 20.4 

+.990 .995 36.69 0.55 53.0 40.67 0.28 36.4 


APPENDIX B - DESCRIPTION OF TOWING TANK AND INSTRUMENTATION 


An early description of the Ship Model Towing Tank, 
prior to the installation of the towing carriage, is given in 
a paper entitled "The Ship Model Towing Tank at M.I.T.,'' Society 
of Naval Architects and Marine Engineers, 1953. This paper 
does, however, give an up to date description of the tank di- 
mensions and the speed measuring and counting equipment. Since 
the towing carriage was installed late in 1961, a brief de- 
scription will be included ‘here. 

The carriage is constructed of aluminum and rides on 
two hollow cylindrical stainless steel rails, one overhead and 
one mounted on the tank's side wall. The overhead rail is 
Simply for support and guidance of the carriage which rides it 
on six fibreglass wheels, two above and one below the rail at 
either end of the carriage frame. The carriage is also sup- 
ported by the side rail with three wheels riding on the after 
end and a soft rubber drive wheel connected via a coupling to 
Preicctric motor on the forward end. The electric drive motor 
is controlled by the console at one end of the tank via cables 
attached to pulleys which slide along a smaller rail mounted 
on the tank wall above the side rail. The control cables are 
free riding and consequently as the carriage moves from one 
end of the tank to the other the speed decreases as more cable 
length is either drawn. out from one end or piled up. A dif- 
ferential pulley system was installed which helped in alleviating 


this situation to some degree. 
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The carriage speed can be controlled either manually by 
setting a dial on the console or automatically at a set speed 
or to position the carriage over a self propelled model. The 
automatic control was inoperative during the period of the ex- 
periments included in this thesis, and speed was controlled 
manually. The speed measuring apparatus includes a type 
500 Hewlett Packard frequency meter and a type 921 Hewlett 
Packard counter described in greater detail in the foremen- 
tioned reference... 

Figure B-I shows a schematic drawing of the model dye 
system which is self-explanatory and figure B-II shows the 


details of the drag measuring towing dynamometer. 
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APPENDIX C - COMPUTER PROGRAMS 


The IBM 7090 electronic computer located at the M.I.T. 
Computation Center was utilized as an aid in reducing drag 
data and in computing the theoretical portions of the thesis. 

Three programs were written: 

Program I = Reduction of drag data. 

Program II- Calculation of turbulent boundary layer charac- 
teristics on the bodies tested. 

Program III- Calculation of the pressure velocity distribu- 
tion on the bodies tested. 

These programs were written in 7090 FORTRAN. On the 
next several pages the following are given for each of the 
three programs: 

1. FORTRAN program listing. 

2. A list of program symbols and their meanings. 

3. A logical sequence flow chart of the program. 

Both FORTRAN and compiled binary machine language card 
decks and listings are available for these programs in the 
Department of Naval Architecture and Marine Engineering 


computer program files. 
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SYMBOL 
RO 

VIS 

SLEN (N) 
WSUR(N , M) 


DT (L,M,N) 


v‘(L,M,N) 
N 


M 
L 
SL 


CFS 


CR 


PROGRAM I - LISTING OF SYMBOLS 


MEANING UNITS 
Mass density of tank water slug/ft? 
Kinematic viscosity of tank water ft2/sec. 
Body length ft. 
Wetted surface of body at 9 
particular depth EC; 
Measured drag lbs. 
Model velocity kts» 
Body. index 

Depth index 

Speed index 

Lenth of strut submerged inches 
Strut drag Lb. 
Corrected body drag Lb. 
Speed - length ratio(V/7 L) writ 
Reynolds number on SLEN (N) 

1/2 §s v? lb. 


Total drag coefficient 


Friction drag coefficient using 
Shultz-Grunow. formulation 


Friction drag coefficient using 
Schoenherr formulation 


Residual drag coefficient 
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LOGICAL SEQUENCE FOR PROGRAM 1 


START 
READ § , (RO) 
and 77 , (VIS) 
READ LENGTH OF 
BODY, (SEN (N)) 


READ MEASURED 
DRAG, DT (L,M,N) 





) 


COMPUTE STRUT DRAG 
(2) (DSTR) , Equation (27) 






DSTR=0.0++JCOMPUTE CORR. 
DRAG, (DTC) 





COMPUTE V/L/2(SLR) 
and REYNOLDS NUMBER, (RN) 







COMPUTE C;, (CT) 
READ VELOCITY on 
Vv (L,M,N) 






i 
\ 








COMPUTE Cf from Schultz- 
Grunow Formula, (CFSG) 
Equation (30) 





SET BODY INDEX 
N=] 


PRINT BODY 
NUMBER = N 


ser DEPTH INDEX 
Mel 












CFSG=CFS 


COMPUTE Cf from Schoenherr 
Formula using CFSG 
Equation (29) 


IF(ABS.VAL. (CFSG-CFS) -.001 


Y 
<G> UF M-2 <=> lcompure Cy, (CR) 
¢ a ae 
STRUT LENGTH, (PRINT V(L,M,N),SLR,DT(L,M,N) ,DTC, 
' |SL=6.0 | CT,CF,CR, and RN 


STRUT LENGTH, 


and OUTPUT HEADING 


@)——= 
Tg DEPTH=M 


STRUT LENGTS ,| 
SL=0.0 





SL=14.0 (IS L=6 L=L+1 
SET VELOCITY INDEX 
ao L=] = > 


zs GF U2 ><> 
f b_}. 
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PROGRAM I FOR CALCULATION OF CR ON BODY 


FORTRAN | 

DIMENSION DT(6,3,5),V(6,3,5) ,SLEN(5) , WSUR(5,3) 

READ 24,RO, VIS 
24 FORMAT (F5.3, E19%.4) 

READ 25, (SLEN(N) ,N=1,5) 
25 FORMAT (5F4.2) 

READ 26, ((WSUR(N,M), N=1,5),M=1,3) 
26 FORMAT (14F5.3) 

READ 27, (((DT(L,M,N), L=1,6) M=1,3),N=1,5) 
27 FORMAT (18F4.2) 

READ 28, (((V(L,M,N), L=1,6), M=1,3),N=1,5) 
28 FORMAT (18F4.2) 

DD 59 N=1,5 

PRINT 32,N 
32 FORMAT (14H1 BODY NUMBER=11) 

DP 50 M=1,3 

PRINT, 33,M 
33 FORMAT (8H@ DEPTH=11) 

PRINT 34 
34 FORMAT (111H6 VEL. (KT) REYNOLDS NO. SPEED/LENGTH TOT.DRA 

1 & TOT. DRAG (CORR) CT CF CR/) 

IF (M-2) 19,11,12 
19 SL=14.9 

GO TO 14 
11 SL=6.@ 

GO TO 14 
12 SL=9.@ 
14 DO56 L=1,6 

IF (M-2) 61,61,60@ 
69 DSTR=9.@ 

GO To 62 
61 DSTR=(RO* ((1.689*V(L,M,N) )**2.0)*((0.964*SL)+1.17))/288.0 
62 DTC=DT(L,M,N)- DSTR 

SLR=V(L,M,N) /SQRTF (SLEN (N) ) 

RN=1 .689*V(L,M,N)*SLEN(N)/VIS 

A=. 5*RO*WSUR(N ,M)* (1.689*V(L,M,N) )#*2.0 

CT=DTC/A | 

X1=9.434294*LOGF (RN) 

CFSG=@.0725/ (X1-2.0)**2.6 
13 X2=9.434294*LOGF (RN*CFSG) 

CFS=9.9585/x2**2.0 

IF (ABSF (CFSG- (¢. @G1LE-93) )7,7,15 
15 CFSG=CFS 

GO TO 13 
7 CR=CT-CFS 

PRINT 35, V(L,M,N),RN,SLR,DT(L,M,N),DTC,CT,CFS,CR 
35 FORMAT (3H F4.2,7H 1PE19.3,9H OPF5.2,1CH 

1@PF5.3,8H @PFS.3 , 16H 1PE19.3,4H  1PE19.3,5H 

21PE19.3) 
59 CONTINUE 

CALL EXIT 

END 
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SLEN (N) 
V(L,M,N) 
X(I,N) 
R(I,N) 


AL(I,N) 


Cr EF 2 


PBLT 


BLTO 


PROGRAM II - LISTING OF SYMBOLS 


Meaning 


Reciprocal of power in 
velocity profile 


Kinematic viscosity of tank water 
Body length 

Model velocity 

Distance from body nose 

Radius of body at X(I,N) 


Cosine of angle between body sur- 
face and body axis at X(I,N) 


Body index 

Depth index 

Speed index 

Location on body index 
Reynolds Rumber on X 


Boundary layer thickness on 


flat plate (d6pl1) 


Displacement thickness on flat 
plate (<* pl) 


Estimate of boundary layer 
thickness on body 


Boundary layer thickness on 
body ( 


Ratio of boundary layer thickness 
to model radius 


An pl 


A\2n pl 


Displacement thickness on body 


(& *) 
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Units 


fto aac. 
£t,. 
kts. 
inches 


inches 


inches 


inches 


inches 


inches 


inches 


' m4 





PROGRAM II (CONT'D) 


Symbol _Meaning Units 
OMNP Ga an 
pl 

OMTNP ir. a -- 
BMT Momentum thickness on body (Q) inches 
RC Body radius corrected to 

account for displacement inches 

thickness 
Y Normalized radius 
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LOGICAL SEQUENCE FOR PROGRAM II 
START 


ee tee a 


READ 1/n, (P);27, (VIS) ;and 
BODY LENGTH, (SLEN(N) ) 





COMPUTE @ ON FLAT PLATE 
(BMT), Equation (3a) 


READ VELOCITY set S ON BODY* SON FLAT PLATE 
V(L,M,N) (BLTO=PBLT ) 


READ LOCATION ON BODY 


COMPUTE 45 ON BODY 
X(1I,N) 


(BLT) ,Equation (1) 


READ BODY RADIUS BLTO=BLT 
R(1I,N) . 


<t>({IF (ABS.VAL. (BLTO-BLT)-.001) 


COMPUTE 4* ON BODY, 
SET BODY INDEX (DT), Equation (2) 
N= | 


(5) COMPUTE 9 ON BODY, 
SET DEPTH INDEX (BMT), Equation (3) 


‘ i] 
Lf M= M+ 
Tad (4) COMPUTE CORRECTED BODY RADIUS, 
= VELOCITY INDEX (RC), Equation (4) 
L=1 






READ COSINE OF ANGLE BETWEEN 
BODY SURFACE AND BODY AXIS 
AL(N, 1) 


LeL+1 


COMPUTE NORMALIZED BODY RADIUS, 
PRINT OUTPUT HEADINGS (Y), Equation (4a) 


PRINT X(I,N),R(1,N),RC,BLT, 
<NoO> DT ,BMT,RN,and Y 
@Q | 


COMPUTE REYNOLD'S IS L=3 «<YES>(IS_1=16) <K>>-(2) 
NUMBER, (RN) 


ComPUTE 5 ON FLAT PLATE,| <YES> 


SET LOCATION INDEX 


‘] mm 





(PBLT) ,Equation (14 ) <Ko> IS M=3 
> 
OMPUTE §* ON FLAT PLATE, OQ 
(PDT), Equation (2b) 





oo 
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PROGRAM II, CALCULATION OF BODY DISPLACEMENT THICKNESS ,ETC. 


FORTRAN 
DIMENSION V(3,3,5),X(16,5),R(16,5),SLEN(5) ,AL(16,5) 
READ 21,P,VIS, (SLEN(N),N=1,5) 
21 FORMAT (F3.1,E19.4,5F4.2) 
READ 22, (((V(L,M,N),L™1,3) ,M=1,3),N=1,5) 
22 FORMAT (18F4.2) 
READ 23, ((X(I,N),1I=1,16) ,N=1,5) 
23 FORMAT (14F5. 2) 
READ 24, ((R(I,N),I=1,16),N=1,5) 
24 FORMAT (16F4.2) 
READ 25, ((AL(I,N),I=1,16),N=1,5) 
25 FORMAT (14F5.3) 
DO 59 N=1,5 
DO 50 M=1,3 
DO 50 L=1,3 
PRINT 1006,N 
199 FORMAT (14H2 BODY NUMBER=T1) 
PRINT 110,M 
116 FORMAT (13H BODY DEPTH=I1) 
PRINT 101,V(L,M,N) 
191 FORMAT (16H0 VELOCITY (KTS)=F4, 2) 
PRINT 102 
192 FORMAT (95H®@ X(INCH) RADIUS CORR. RADIUS DELTA’ DELTA 
4STAR THETA REYNOLDS NO, Y-NORM/) 
DO 50 I=1,16 : 
RN=1.689*V(L,M,N)*X(1,N)/(VIS*12.9) 
PBLT=@ .38*X(1I,N) /RN**Q.2 
PDT=PBLT/(1.9+P) 
BLTO=PBLT 
12 BLT=PBLT/(1. 6+ (BLTO/(3.0#R(I,N))) )##9.8 
IF (ABSF(BLTO-BLT)=-.601) 7,7,1¢ 
16 BLTO=BLT 
GO TO 12 
7 A=BLT/R(I,N) 
DLNP=PDT/PBLT 
DLTNP=DLNP* (P+1.6)/ (2. G*P+1.0) 
DT=PDT* (1. 6+ A-SQRTF (1 , §+2. G#A* (1. G—DLNP)+ (A**2.9)*(1.0-DLTNP)))/(A 
1 *DLNP* (1. 0+A/3.0)**6.8) 
OMNP=p/((P+1.@)* (P+2.9)) 
OMTNP=OMNP* (P+2.0)/2.0*P+1.@) 
BMT=BLT* (OMNP+A*OMTNP/2.0) 


RC=R(1I,N)+AL(1,N)*DT 
Y=RC/ (SLEN (N) *6 .@) 
PRINT 193,X(1,N),R(1,N),RC,BLT,DT,BMT,RN,Y 
193 FORMAT (3H F5,2,6H F4.2,7H _ F5.3,9™ F5.4,7H 
1 F5.4,6H F5.4,5H 1PE19.3,4" OPE12.5) 
50 CONTINUE 
CALL EXIT 
END 
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Symbol 
Bo 
Y(1) 
R(1) 

RO 

PATD (M) 
V(L,M,N) 
N 

M 

L 

AL (1) 


A(1) 
B(I) 
G(J,1) 
C(J,1) 
D(J ,0) 
SUM(I) 
E(K,1) 
H(TI) 


PROGRAM III - LIST OF SYMBOLS 


Meaning Units 


Gaussian abscissae, number (1-16) 

Normalized Gaussian abscissae based on model 
length of two, from -1 to+ 1 

Normalized ordinate of body meridian profile 


Gaussian weighting functions 


Water density slug/ ft3 
Static pressure (gage) PsiG 
Absolute model velocity kts. 


Body number 
Depth 
Speed 


Cosine of angle between tangent to body profile 
and x-axis 


A(1) ratio of length to diameter 

k, (1), longitudinal virtual mass coefficient 
g(x,t), equation of prolate spheroid whose ends 
coincide with ends of the given body, and which 
intersects body at x =t. 

K (x,t), matrix of functions of body 

K’ (x,t) ,Matrix of functions of prolate spheroid 
Z(K=K*) 


E(n,t),error functions 


kK, + E 
1 +k (x) 


Iteration number 


G2 





PROGRAM III = (CONT*D) 


Symbol Meaning | Units 

SUM (I) 2 (Bn, - Ey) 

0(1) =] Eni 

U(T) Va ; ratio of flow velocity over body to free 
stream velocity. 

P(T) p/a , pressure distribution 

VEL(I) Absolute flow velocity over body kts. 

PCC (TI) _Gage pressure on body lb/ft? 
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LOGICAL SEQUENCE FOR PROGRAM III 


Read X(I), R(I), zO, 
PATD(M), V(L,M,N) 


t= N,M,L, AL(17) 








Read Y(1),AL(I) 





Complnte H(I) = 
(8) B(17) + EQ,1) 
iL eG 










Compute ZK(E_.- E 
(EK (J, 1)x(E (Kd) -E(K 












Compute A(x); (A(T)) 
(4 » equation (17) 







equation 
Compute B(17) 








Compute g(x,t), | Compute E,.) 4 P =P +1 
(G(J,1)) eq. (16) (E(K+1,1T)) 
eq. (20) (p ) 
Compute K(x,t), 
(C(J,1) eq, (14) I= I +1 
Compute K* (x,t), | 9A 
(D(J,1)) eq. (15) 


(comtinued on next page) 


G4 
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Compute U(I) 
=(1+k) +2E,;) cos 4 


= (14B(17)+ZE(K,1I))AL(I) 


Print N,M,L and 
output heading 


Compute E,j; 
E(K,I),K= 1..5 





















Compute py/q 
P(1)=1_U(1)2 













Compute 


VEL (1)=U(19EV(L, M,N) 


> m2 


Compute Poo. 
PCC (I)=P (1) xROXV (L ,M,N)2 













re. 
+ PATD(M) 


[Print X(I),Y¥(I),UCI), 
iP (I) ,PCC (1) , VEL(T) 


Is N= 5 


, 2 


"Sa tes 


yes 
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PROGRAM III CALCULATION OF PRESSURE DISTRIBUTION 


FORTRAN 
DIMENSION X(16),R( 16) ,Y(16) ,AL(16,5) ,A(17) ,B(17) ,C(16,16) ,D(16,16) 
1,SUM(16) ,E(5, 16), H(16),0(16) ,U(16), $(16) 'P(16), G6, iia), PATD(3), 
ev (3, 3,5), VEL (16), PCC (16) 
READ i) (x(Z), R(1), T=1,16) 
1 FORMAT (6E12.5) 
READ2,, RO ,;’(PATD(M) ,M=1,3) 
2 FORMAT (F5.3,3E11.4) 
READS, (((V(L,M,N) ,L=13) ,M=1,3),N=1,5) 
3 FORMAT (18F4. 2) 
READ7, ((AL(I,N) ,I=1,16) ,N=1,5) 
7 FORMAT (14F5.3) 
4 READS,N,bi,L,A(17) 
5 FORMAT (311. E12. 5): 
READG, (Y(I) ,I=1,16) 
6 FORMAT (6E12. 5) 
DO 19 I=1,16 
A(T) =SQRTF ((X(T)41.)*(1. ~X(I))/(¥ (1) **2. )) 
19 CONTINUE 
D029 r=1,17 
WD=(A(I)*LOGF(A(1)+SQRTF(A(I)** 2 »-1.) )-SQRTF(A(I)**2.-1.))/. 
1( (A(T) **2. )*SQRTF (A(T) **2.=1. )-A(L)*#LOGF (A'(L)+SQRTF(A (I) **2., -1.))) 
20 CONTINUE 
DO 30 I=1, pi 
DO 36 J=1 
G(J, re (Liz) *ee, )* (X(J)41.)*(1.-X(J) )/ ((X (1) 41.) * (1. -X(T))) | 
C(J, “T)=R(S)* (¥(S)**2. Y/ CC CX GI) —X (I) )¥#2 .4Y (J) #*2. )#*1, 5) 
D(J ,T)=R(J)*G (J ,T)/(((X (JS) =X (1) )#*2.4G(S , 1) )**1.5) 
39 CONTINUE 
DO 4¢ I=1,16 
SUM(I)=@. G 
DO 39 J=1,16 
S(J)=C (J, T)-D(J, T) 
SUM(I)=SUM(I)+S (J) 
39 CONTINUE 
K=] 
E(K,I)=1.-(((4.4B(17) )/2. )*SUM(I) (1.4B(17))/(1.+B(T)) 
H(T)=(E(K,1)+B(17))/(1.+B(17)) 
40 <i ln 
DO 5@ K=1,4 
DO 5¢ I=1,16 
SUM(1)=@. " 
DO 49 J=1,16 
S(J)=C(J, T)*(E(K, J)-E(K,1)) 
SUM(I)= SUM (1) 4 S(J) 
49 CONTINUE 
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59 
52 
53 
54 
55 


PROGRAM III (CONT'D) 


BE (KHT)=H (I) *E(K,I)-.5*SUM(T) 
CONTINUE 

PRINT 52,N 

FORMAT (14H2 BODY NUMBER=-T1) 

PRINT53 , M 

FORMAT (8H% DEPTH=11) 

PRINT54 ,V(L,M,N) 

FORMAT (8HG SPEED=F4., 2) 

PRINT55 

FORMAT (83H@ X(T) Y¥(T) U(T) P(T) 
1 PCC (TL) VEL (1) /) 

DO 66 I=1,16 

O(T)=E(1 , 1)+E(2,1)+E(3,1)+E(4,1)+E(5,1) 
U(I)=(1.4B(17)+0(1) )*AL(T,N) 

P(I)=1.-(U(I))**2. 

VEL(I)=U(I)*V(L,M,N) ~ 

PCC (I) =P (1)* (. 5*RO* ((V(L,M,N)*1.689)**2. ) )4+PATD(M)*144., 
PRINT 59, (X(I),Y(1),U(1), PCI) ,PCC(1), VEL(Z)) 

FORMAT (6E14.5) 

CONTINUE 

GO TO (4,4,4,4,61),N 

GO TO (4,4,62) ,M 


GO TO (4,4,63) ,L 
CALL EXIT 
END 
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APPENDIX D — EXPERIMENTAL DATA 


TABLE D-I ~ EXPERIMENTAL DRAG DATA 


Dial 


20 
20 
24 
24 
28 
28 
20 
20 
24 
24 
28 
28 
20 
20 
24 
24 
28 
28 
20 
20 
24 
24 
28 
28 
20 
20 
24 
24 
28 
28 
20 
20 
24 
24 
28 
28 


Speed Voltmeter Voltmeter Voltmeter Measured 


Run Stern Depth Speed (KTS) Zero(V) AVG (V)  Diff.(V) Drag(lbs) 


0.80 
0.77 
1.12 
IztZ 
1.48 
1.48 
0.75 
0.74 
1.09 
1.09 
1.50 
1.51 
0.77 
0.79 
1.16 
1.17 
1.62 
1.58 
0.77 
0.77 
1.16 
1.21 
1.58 
1.58 
0.75 
0.75 
1.19 
1.11 
1.61 
1.61 
0.74 
0.74 
1.13 
1.16 
1.56 
1,53 


-43.7 
-45.0 
~46.5 
~45.7 
-45.0 
-43.5 
-45.1 
~47 .7 
~45/5 
-44.7 
-43.0 
~43.5 
-46.0 
~46.1 
-45.5 
~45.0 
~46.0 
~45.0 
-45.0 
~46.0 
~46.0 
-46.5 
-46.0 
~45.5 
-35.0 
-36.0 
-35.0 
-37.5 
-42.5 
~44.0 
-25.0 
~-26.7 
-39.0 
-40.5 
-44,0 
-46.0 
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-~26.4 
-27.9 
-20.4 
-18.5 


~14, 75 


-11.7 
-29.2 
-30.8 
~22.6 
-18.5 
-18.4 
-19.0 
=30.7 
-30.9 
-23.6 
~25.1 
~19.6 
-21.0 
-23.2 
-21.2 
-24.0 
-25.8 
-19.7 
-20.0 
-22.6 
-24.2 
-13.2 
-~15.2 
-12.5 
-12.6 
~15.4 
-16.1 
-16.1 
-15.2 
=09.2 
~12.6 


17.3 
17.1 
26.1 
27.2 
30.25 
31.8 
15.9 
16.9 
22.9 
26.2 
24 .6 
24.5 
15.3 
15.2 
21.9 
19.9 
26.4 
24.0 
21.8 
25.0 
22.0 
20.7 
26.3 
25.9 
12.4 
11.8 
21.8 
22.3 
30.0 
31.4 
09.6 
10.6 
22.8 
25.3 
34.8 
33,4 


0.21 
0.21 
0.33 
0.34 
0.38 
0.40 
0.20 
0.21 
0.24 
0.33 
0.31 
0.31 
On19 
0.19 
0.27 
0.25 
0.33 
0.30 
0.28 
0.31 
0.28 
0.26 
0.33 
0.32 
0.16 
0.15 
0.28 
0.28 
0.38 
0.39 
0.12 
0.13 
0,28 
0.31 
0.44 
0.42 





TABLE D-I (CON’T) EXPERIMENTAL DATA 


Dial Speed Voltmeter Voltmeter Voltmeter Measured 
Run Stern Depth Speed (KTS) Zero(V) ‘AVG (V) Diff.(V) Drag(1bs) 


mo 2 20 0.75 -44.0 -28.6 15.4 0.19 
19A 20 0.73 -45.0 -29.4 15.6 0.20 
20 24 1.14 -46.0 -20.0 26.0 0.32 
20A 24 1.16 -46.5 -17.9 28 .6 0.36 
21 28 1.61 -44.5 see 31.3 0.39 
21A 28 1.56 -43.5 ~17.9 25.6 0.32 
0.80 -24.0 -09.3 140a 0.16 

2 0.73 -24.5 =1103 13.2 0.14 
1.15 -28.0 -09.0 19.0 0421 

1.14 -26.0 -07.5 18.5 0.21 

1.67 -25.0 -04.6 20.4 0.23 

1.60 -24,5 -04.0 20.5 0.23 

0.43 -37..0 mien 23.8 0.28 

0.75 -32.0 -10.3 PPT 0.26 

1,18 —32,0 -09.4 22.6 0.26 

1.18, -30.0 -07.6 22.4 0.26 

1.53 -29.0 -05.5 23.5 0.27 

2 1.50 -26,0 -03.8 22.8 0.26 

4 0.70 -40.0 -08.6 21.4 0.25 
0.70 -42.5 ies 25.0 0.29 

1.00 -45.0 -19,2 25.8 0.31 

1.25 -42.5 -16.2 26.3 0.31 

1.65 -40.0 Sli 28.8 0.35 

1.67 -40.0 alle 28.8 0.35 

0.78 -35.0 Sl7ee 17.8 0.20 

0.78 -34.5 Stone 18.8 0.21 

4 1.22 +31.5 -08.8 2207 0.26 
1.23 -32.0 -10.0 22.0 0.25 

1,63. =S250 -07.4 24.6 0.29 

1.58 -30.0 -05.2 24.8 0.30 

0.77 -43.0 =19 57 2one 0.26- 

0.76 -48.0 — 2A 1 23.9 0.28 

1,20 —44/,0 -15.2 31.8 0.39 

1.18 -40.5 123 28.2 0.34 

1.58 -41.0 -05.4 35.6 0.45 

4 1.58 -44.0 -07.7 36.3 0.45 

5 0.78 -37.0 18) 1 18.9 0.21 
0.80 -44.0 -20.9 23.1 0. 27 

| 1.21 -30.5 -05.0 25.5 0.30 
1.18 -34.5 -06.2 28.3 0.34 





TABLE D-~{(CONT'D) EXPERIMENTAL DATA 





Dial Speed Voltmeter Voltmeter Voltmeter Measured 
Run Stern Depth Speed (KTS) Zero (V) AVG (V) Diff.(V) Drag(1bs) 
39 | 28 1.67 -37.5 -@1.1 36.4 0.46 
39A 1 28 1.80 -43.0 -07,3 35.74 0,44 
40 2 20 0,80 -45.0 -20,1 24.9 0.29 
40A 5 20 0.80 -~46.0 ~20.3 29.4 0.30 
41 24 1.28 ~46.5 -~15.3 oy a 0.38 
41A 24 1.28  ~-46.0 ies 32.5 0.40 
42 28 1.70 -~46.0 ~ 8.8 Slez 0,47 
42A 2 28 1.65 -~46.0 ~ 8.4 37.6 0.47 
43 3 20 0.75 ~43.5 -~18.1 25.4 0.30 
43A 20 0.75 ~-39.5 ~17.3 22.2 0.26 
44 24 1.25 -39.0 -14,4 24.6 0.29 
44A 24 1,28 -38.5 -13.1 25,4 0.30 
45 28 1.63 -37.0 -~07 .4 29.6 0,35 
45A 5 3 28 1.60 ~36.0 ~07.0 29.0 0.35 


Depth 1 = 18" to Keel (4.5*Dia.) 
Depth 2 = 10" to Keel (2.5*Dia. ) 
Depth 3 = 3.6"to Keel (0.9%Dia.) 


Runs - 1 = 21A Conducted 24 March 1962 
22 =~ 45A Conducted 25 March 1962 


Water Temp = 51°F 


rs) 


U = 1.386 x 10° ft/sec 


| 4 
S = 1,939 lb - sec“/ft 


Dial speed is setting on manual speed control knob. 
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TABLE D-II —- SEPARATION POINT DATA 


Run Stern Run Separation Body 
Classification Point Angle 
1 3 LD 5.0" 33° 
1A LD 5.0 33 
2 MD 5.0 33 
2A MD 2.0 33 
3 HD 5.0 33 
3A HD 5.0 33 
4 LM 4.75 28 
4A LM 4.75 28 
5 MM 4.75 28 
5A MM 4.75 28 
6 HM 4.75 28 
6A HM oe 33 
7 LS 4.5 23 
TA LS 4.5 23 
8 MS 4.5 23 
8A MS 4.5 23 
9 HS 4.75 28 
9A 3 HS 4.5 230 
10 1 LS 0.8 24 
10 LS 0.8 24 
11A MS 1.25 A2 
11B MS 1.25 A2 
12 HS 1.3 42 
12A HS g eS 5 42 
13 LM 0.75 24 
13A LM 0.75 24 
14 MM 1.0 32 
14A MM ow SZ 
LS HM bel 32 
15A HM 1.1-1.2 32 
16 LD 0.8—09 24 
16A LD 0.8 24 
17 MY) 0.8 24 
17A MD C.8-0.9 24 
18 HD is 42 
18A 1 HD Des 42 
19 2 LD 2.75 30 
19A LD 250 25 
20 MD Jee 36 
20A | MD Sel 36 
21 HD ? = 
21A HD 3.0 36 
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Stern 
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TABLE D-II (cent'd) 


Run 
Classification 


LM 
LM 
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Separation 
Point 


>= 


o © @ @ 


on 


° 


COON OANNNNRKFRFWOWRFFOON WANK OO”) DOO © 


© © e © ° Qe ° e © © © ® e © 


° ° e © @ © 


VE VAMOMMNAVNYVITIYVYYVNVANYIYNWWWWNNNHVNNNW 


© 
2) 
~~ 


Body 
Angle 


36° 
25 
30 
30 


30 
33 
30 
36 
42 
42 
42 
38 
38 
44 
44 
49 
49 
44 
44 
49 
49 
49 
49 
28 
28 
28 


~30 


38 
23 


23 
23 
23 
23 
30 
30 


30 


TABLE D-II (Cent ‘d) 
NOTE: 
Run classification code ~ lst letter designates speed 


L = low 
M= medium 
H = high 


2nd letter designates depth 
D = depth 1 (deep) 
M= depth 2 (medium) 
S = depth 3 (surfaced) 


Separation point is measured in inches from end of cylindrical 
section. Marks on photographs are 4 inch apart. 


Body angle is angle between body axis and body surface at 
separation point. 
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APPENDIX E - DISCUSSION OF DRAG FORCE MEASUREMENTS 


The experimental data from Table D-I was reduced to co- 
efficient form using the following procedure. 
(Drag measured) - (Strut drag) = (Total corrected drag) 
where strut drag was calculated by the method that follows. 
Measurement of a strut tare drag was attempted but due to 
either the end effects on the strut because of the lack of 
a dunmy body or faulty instrumentation these values were 
of a larger magnitude than the values of measured drag of 
the body plus strut and thus were not valid. An engineer- 
ing value of strut drag composed of profile drag, spray 
drag, wave drag, and interference drag with the body was 
calculated as follows: 


Profile Drag: 
Me es 8 = .246 for strut used. 


Reynold Number on Strut ¥* 3 x 104 
from figure 126, page 247 of (9) 
Cd = 0.08 


Drag = q Cd(.8 x SL) = q (.064) SL 
144 144 


where SL = strut length submerged (inches) 


eet 
a7 10 Vy (1.689) ° 


Spray Drag: 1 


Froude number on strut * 4 x 10° 


t = .8 = .246 


from figure 25, page 10-13 of (10) 
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Ca, = .24 (t)"=.088 x = 1,625 
x 


Drag = q Cd, (x) “4: (.1535) 


Wave Drag: 


from figure 24, page 10-13 of (10) 


Cd, = 1.4 


Drag = q Cd, 1? =q (.895) 
4 


144 


Interaction Drag: 


from figure 23, page 8-10 of (10) 


Cot = .t5 x t/c 
Drag = Ca, q t? «= q (,121) 
144 


ES, 


Total Strut Drag: 


2 
Strut Drag = 1 pV, (1.689) | 06481 + 1.17] (27) 
2 144 


The total corrected drag was converted to coefficient 
form by the following relationship: 
C+ = Drag (28 ) 
1/2 ° YK (1.689)2 (5S) 
where S = wetted surface of the model. 
The coefficient of frictional resistance, Cr, was nex calcu- 
lated using Schoenherr’s relationship that: 
Cp= 0.0585 or (R,C4)|-? 
= 0.0585 (29) 


2 
[-4343 loge (Rn Cs) | 
In order to use the electronic computer to evaluate (29), 


the Schultz—-Grunow formulation 
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_2 
Ce = 0.0725 |1ogy, Rn -2.0] 


SS0n07Zo 
5 (30) 
[4343 logeR, -2.0 | 


which yields Cf values greater than that of (29) in the range 
of interest in this problem was calculated and used as a first 
value to solve (29) which is implicit in Cf. 

The residual resistance coefficient, Ci. was then 
calculated from: 


oe = Cc. - C (31) 


f 


It was hoped that perhaps some correlation between the 
change of location of separation on the body and change in C, 
for various runs could be found if our drag measurements were 


good. 
An IBM 7090 commuter program, Program I was written to 


solve equations (27) through (31). The details of this program 

are given in Appendix C. Sample outputs of this program are 
included in Table E-~I at the end of this appendix. These values 
which are typical of the values obtained on all the bodies are much 
too large for the underwater shapes tested. At submerged depth 

1 the values of Cr and Cy are over ten times as large as would 

be expected by comparison to Series 58 data. 

At first strut drag was suspected but in the surfaced con- 
dition where strut drag does not enter into the calculation the 
residual drag values are much too large also. The bow waves 
generated in the surfaced condition are shown in Figure Ix. It 


is not possible for the wave making resistance created in this 
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condition to account for such large values of Cc... Besides 
being too large in magnitude the values are not consistent 
from body to body with the highly streamlined body some- 
times showing a higher C.. value than a less streamlined body 
for similar speed and depth of submergence. Plots of C.. 
versus both Reynolds number and speed-length ratio were 
prepared but these were completely inconsistent and unsatis- 
factory when compared to similar plots of similar bodies in 
the literature. 

At the time of this writing the cause of the excess 
drag is unknown. It is speculated that the dynamic response 
of the drag measuring system is causing a false reading of 
drag due to some phenomena that does not appear during sta 
tic calibration. The calibration curves used in this oa tic 
are shown in Figure E-I. It is also of interest that a= 
nother group using another Shoe model experienced the 
same difficulty, i.e., too much drag. 

If it is desired to make further drag measurements 
on the model it is recommended that the drag force dyna- 
mometer be mounted at the model rather than in the strut. 
This will eliminate the problem relating to strut drag with 
the exception of the interference drag between the strut and 
body. It is also recommended that the dynamic response of 
the entire drag measuring system used during the tests be 


given a thorough investigation. The authors realize that 
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the new towing tank and dynamometer equipment was put to 
use on an emergency basis, the construction of such appara- 
tus being completed, and before the carriage and instrumen- 
tation were properly calibrated, accepted, and put into 
regular use. Upon completion of all of these steps it 
should be feasible to carry out successful drag measure- 


ments on submerged bodies at the new M.I.T. Towing Tank. 
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Figure ErI- Dynamometer Calibration 
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TABLE E-I — SAMPLES OF COMPUTER OUTPUT, PROGRAM I 


BCCY NUMEER=1 


VEPTE=1 
VEL.(KT. REYNOLDS NO. SPEED/LENGTH TOT.DRAG TOT.DRAG (CORR) 
0.74 2.372E U5 Cua 0.120 C.098 
C.74 2.372t C5 C.46 0.130 0.108 
1.13 3.622 U5 Geno 0.280 0.229 
1.16 3.718¢ 05 Ose 0.310 0.257 
1.56 5.000t U5 0.96 0.440 0.343 
1.53 4.904E U5 C.94 0.420 0.327 
DEPTR=2 
VEL.(KT. REYNOLDS NC, SPEED/LENGTH TOT.DRAG TOT -ORAG (CORR) 
0-75 2-404E 65 C246 C. 160 0.143 
0.75 2-4C4E 65 0.46 0.150 0.133 
Lely Ze814E C5 0.73 0.280 0.238 
Bett 3.558b OS 0.68 0.280 0.243 
1.61 5e160E U5 C.99 0.38C 0.363 
L.61 5.160t 05 0.99 0.390 0.313 
DEPTR=3 
VEL.(KT. REYNULUS NO. SPEEO/LENGTH TOT. DRAG TOT .ORAG (CORR) 
0.77 22468t US 0.47 0.280 0.280 
GavT 2.468E 65 Gea? 0.310 Os310 
Le16 3.718E O5 Gave 0.280 0.280 
1.21 3e@7BE C5 C.75 0.260 0.260 
1.58 5.C64t US C.97 0.330 0.330 
1.58 5.C64t U5 G.97 0.320 0.320 
CT CF CR 
2-541E-02 5. 909E-03 1.950E-02 
2.799E-02 5.909E-03 2.208E-02 
DEPTH = 1 2.543E€-02 5 400E-03 2.003E-02 
2. 700E-02 5.371E-03 2.163E-02 
1-998E-02 5 .053E-03 1 -493E-02 
1-978E-02 5 0736-03 1.471E-02 
CT CF CR 
3.604E-02 5.892E-03 3.015E-02 
3.353E-02 5.892E-03 2.763E-02 
DEPTH = 2 = 2.377£-02 5. 342E-03 1.842E-02 
2.795E-02 5 -421E-C3 22253E-02 
1.653E-02 5 -021LE-03 L.-L51LE-02 
1. 707E-02 5 «021E-03 1.205E-02 
oR CF CR 
7.797E-02 5. 858E-03 7-21 1E-02 
8.632E-02 5.85 8E-03 8.047E-02 
DERTH = 3-336 436E-02 5371E-03 2-898E-02 
2.932E-02 56 324E-03 2-400E-02 
2-182E-02 5 .040E-03 1 678E-02 
2.116E-02 5.040E-03 1.612E-02 
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BODY NUMBER=3 


CEPTR=1 
VEL.(KT. REYNGLDS NO. 
0.86 2.856E U5 
On? 2-749E 05 
Pelz 3.999E G5 
be 2 3.999E 65 
1-48 5.284E US 
1.48 5.284t 05 
DEPTh=2 
VEL. (KT. REYNOLDS NO. 
0.75 2.678E O05 
OA 2.642t 05 
1.09 3.892t U5 
1.09 3.892E U5 
le oc 5.356t US 
151 5.392E 65 
DEPTH=3 
VEL.(KT. REYNOLUS NO. 
a7? 26749t US 
O79 2.821t U5 
1.16 4e142E 65 
1.17 4.178 65 
1.62 5.784t 05 
1.58 S.641E 05 
CT 
3-771E-02 
4.111E-02 
—DEPTH = 1 2.920E-02 
3.024E-02 
1.749E-02 
1.868E-02 
CT 
4.258E-02 
4.623E-02 
DEPTH = 2 2.801E-02 
3.241E-02 
1.411E-02 
1.387E-02 
CT 
5.058E-02 
4.805E-02 
DEPTH =3 3.167E-02 
2.882E-02 
1.985E-02 
1.897E-02 


TABLE E-I_ (CON'T) 


SPEED/LENGTH 


C.47 
C.45 
C.65 
C.65 
C.86 
C.86 


SPEED/LENGTH 


C244 
C.43 
0.64% 
C.64 
0.88 
C.88 


SPEtD/LENGTH 


C.45 
C.46 
C.68 
0.68 
C.95 
C.92 


CF 
5-678E-03 
56 724E-03 
5«290E-03 
56 290E=03 
4.997E-03 
4.997E-03 


CF 
Ds (966703 
Se (T3E-03 
5.320E-03 
5.320E-03 
4.983E-03 
4.976E-03 


CF 
De 724E-03 
5.693E-02 
Sec o2E—O03 
332436203 
4.906E-03 
4.931£-03 
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TOT. ORAG 


G-210 
0.210 
0.330 
0.340 
0.38C 
0.400 


TOT.ORAG 


0.200 
0.210 
0.290 
0.330 
C.310 
0.310 


TOT.DRAG 


0.190 
0.190 
0.270 
0.250 
0.330 
0.300 


CR 
36203 -—02 
53653 9E=02 
26391E=-02 
2-495E-02 
1.249E-02 
1.369E-02 


CR 
3.682E-02 
4.046E-02 
2-269E-02 
2-709E-02 
9-126E-03 
8.894E-03 


CR 
4.485E-02 
4e232E—02 
2-642E-02 
Zed 20e-O2 
1.494E-02 
1.403E-02 


TOT.DRAG (CORR) 


0.185 
0.186 
0.280 
0.290 
0.293 
0.313 


TOT.DRAG (CORR) 


C.183 
0.194 
0.255 
0.295 
0.243 
0.242 


TOT~DRAG(CCRR) 


0.190 
0.190 
0.270 
0.250 
0.330 
0.300 





BODY NUMBER=5 


CEPTR=1 
VEL.(KT. REYNOLDS NC. 
0.78 3-241leE G5 
0.8 3-324 C5 
1.21 DeO26E O05 
1.18 4.903E O5 
1-67 6.-94Gt U5 
1.80 7-480E U5 
DEPTH=2 
VEL.(KT. REYNCLUS NO. 
0.8¢C 34324F G5 
0.8C 32324t U5 
1.28 5e319F C05 
1.26 Se 319E G5 
1.70 71.064E- 05 
1.65 6.-857t O05 
DEPTH=3 
VEL.(KT. REYNCLOS NQ. 
0.75 Z3elk7t O5 
0.75 @eil7tE O05 
1.25 5e194t G5 
1.2% 52319t O5 
1.63 6.773 O05 
1.60 6-649t O5 
CT 
3.540E-02 
4-429E-02 
DEPTH = 1 1e-915E-02 
2-370E-02 
1.452E-02 
1-114E-02 
CT 
-4.905E-O02 
5.-086E-02 
-DEPTH = 2 2-342E-02 
2-483E-02 
1.539E-02 
1.655E-02 
CT 
72374E-02 
6-391E-02 
DEPTH = 3 2.566E-02 
2-532E-02 
1.821E-02 
1.890E-02 


TABLE E-I_ (CON'T.) 


SPEED/LENGTH 


0-42 
0.43 
0.66 
C.64 
C.90 
C.97 


SPEED/LENGTH 


Q.43 
0.43 
0.69 
C.69 
G.92 
C.89 


SPEED/LENGTH 


0.41 
0.41 
C.68 
0.69 
C.88 
0.87 


CE 
5-528E-03 
5+498E-03 
5-04/7t-03 
5-073E-03 
4.731E€-03 
4.662E-03 


CF 
5. 498E-03 
5-498E-03 
4 990E-03 
4.990E-03 
4.714€-03 
4-742E-03 


CF 
5-574E-03 
5.014&-03 
4 T11E-0: 
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TOT .DRAG 


0-210 
0.270 
0.300 
0.340 
0.460 
0.440 


TOT .DRAG 


0.290 
0.300 
0.380 
0.400 
0-470 
0.470 


TOT -ORAG 


0.300 
0.260 
0-290 
C.300 
0.350 
0.350 


CR 


2-987E-O02 
3-879E-02 
1.863€-02 
9.785E-03 
6-478E-03 


CR 
4.355€-02 
4.537E-02 
1.843E€-02 
1-984E-02 
1.067E-02 
1.-181E-02 


CR 
6-81 7E-02 
5 -834E-02 
2-065E-02 
2.033E-02 
1.346€-02 
1.413E€-02 


TOT.DORAG (CORR) 


0.186 
0.245 
0-242 
0-285 
0.349 
0.311 


TOT.DRAC(CORR) 


0.271 
0.281 
0.331 
G.3aor 
0.384 
0.389 


TOT -DRAG (CORR) 


0.300 
0.260 
0.290 
0.3C0 
0.350 
0.350 


APPENDIX F = SAMPLE CALCULATIONS 
Calculation of Hydrostatic Pressure. 


The hydrostatic pressure to the keel of the bodies at 
the three depths of submergence is desired for use in com- 


puter Program III to compute the gage pressure on the body 


from the following relationship: 


4 


where P = dynamic pressure (ps£) 
q=1/20U ? Ost ) 


Ps = hydrostatic pressure (psi) 


Pgage = (2)3 + Ps (144) 


Ps = SJh 
172 
where 9° = 1.935 lb.sec? at S1°F 
ft 
Q = 32.174 ft/sec” 
nN = depth to keel (inches) 


Using the foregoing relation the following results were 


obtained: 
Depth h (inches) Ps (psi) 
1 * - 18 .6484 
2 10>, . 3602 
3 3.6 1297 
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Calculation of Point of Separation from Stratford Criterion 


A typical calculation of the location of the separation 
point according to Stratford's criterionis given below for body 
#3 at velocity of 0.75 knots. 

From the form of (21) separation can occur only when 
Cp , as plotted in Figure F-II, is positive since negative 


values of Cp will result in imaginary numbers. 
At X = 34.8" X= 34.8 =0.99 
L 
from Figure F-I; U = 0.75 
Uoo 
from Table G-I; Reynolds number = 2.65 x 10° at X=34.8" 


thus Rs = .75 (2.65 x 10°) = 1.99 x 10° 


from Figure F~II; Cp = 0.05 and dCp = 0.50 
dx 


n = log 10 Rs = 5.30 from (22) 


since, a“p 2 0; B= .73 from (24) 
dx? i" 
thus the left hand side of (21) becomes 


5.30=2 PM 
(2x.05) 4 (34.8 x .50)° = .625 


while the right hand side of (21) becomes 


-5 _ 6.71/10 
(1.06) (.73)(1.99 x 10°°x10") = 912. Thus the 


criterionis not satisfied at X = 34.8" since the RHS is much 
larger than LHS. Next try 
“ 
X= 34.9 K= .995 
L 
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from Figure F-II; Cp = 0.10 and dCp = 0.96 
dx 


thus the LHS of (21) becomes 


3 


(2 x16) (34.9x.96)” = 1.46 which now makes the 
RHS smaller than LHS so the point of separation lies between 
X= 34.8" and X = 34.9." 

Calculations similar to the above ae used to cal- 
culate the separation points on the other bodies at the va- 
rious speeds. During these calculations it was noted tha the 
criteria was insensitive to speed changes and that the same 
point of separation occurred for all values of speed. The 
points of separation according to this criterionare marked 


by small circles on Figure F-II. 
Calculation of Point of Separation from Buri’s Criterion 


A typical calculation of the location of the separation 
point according to Buri'’s criterionis given below for body #3 
at a velocity of 1.12 knots, 


From the form of equation (26) separation will occur 


only when df U\ as determined from Figure F-I is negative. 
Uoo 


a = -—,115 and U = 1.103; from Figure F-I 
__\Uoo ” Uoco 


and 9 = .091; from Table G-I 
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thus the left hand side of (26) becomes 

-(,115) (1.103) (.906) (091) = -- .0105 7 ~—,006 
Thus the point of separation lies slightly forward of x = .906 
but in the region of where 

"(aos) is still negative. Separation point is at 

Uoo 
X= .900 or X = 31.7 inches which is as close as 
the point can be calculated within the accuracy of the plot in 
Figure F-I. 

Calculations similar to the above were used to calculate 
the separation points on other bodies at various speeds. During 
these calculations it was noticed that the criterion was insensitive 
to speed changes and that the same point of separation occurred 
for all values of speed. The points of separation according to 


this criterionare marked by small circles on Figure F-lI. 


Figures F-I and F-II on the next two pages were plotted 
from values obtained from the outputs of Programs II and III 


as aids to evaluate (21) and (26). 
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APPENDIX G = SAMPLES OF COMPUTER OUTPUT, PROGRAMS fland ITT 


BOOCY NUMRER=3 


BODY DEPTH=1 


VELOCITY(KTS)=1.12 


XCINCH) 


0.19 


wO.o7, . 


2-36 
4.30 
6.72 
9.53 
12.62 
15.90 
19.24 
22a 32 
25262. 
26.43 
30-85 
32.79 
34.17 
34.96 


RADUIS 
0.40 
0.94 
1.37 
1.67 
1.87 
1.97 
2.40 
2-00 
2-00 
2-00 
2200 
2-00 
1.94 
1.6? 
1.15 
Ued5T7 


KQDY NUMBER =3 


BORY -LEPTR=2 


VELCCITY(KIS)=1.51 


XC INCH) 


0.19 
0297 
2.36 


430° 


oan a4 

9.53 
12.62 
15.90 
19.24 
22.52 
25.62 
28.43 
50.85 
32.79 
34.17 
34.9 


KADUIS 


0.40 
0.94 
ba3T 
1.67 
1.37 
laa? 
2.00 
2.00 
2.00 
2-00 
2-00 
2.00 
1.94 
1-62 
Ll 
0257 


CORR. 


TABLE GeI 


0.401 
0.946 
1.383 
1-691 
1.900 
2.003 
2-047 
2-054 
2.061 
2.068 
2-073 
2.078 
Zate | 
1.697 
bea2iy 
0.605 


CORR. 


0.401 
0.946 
1.383 
1-690 
1-848 
2-006 
2-044 
2-052 
2.058 
2-065 
2-070 
Zen > 
2.018 
1.694 
1.214 
0.604 


RADIUS 


RADIUS 


DELTA 
20154 
0564 
21141 
- aa? 
62599 
- 3408 
04225 
25032_ 
5806 
26528 
-7181 
-7751 
8206 
8427 
8326 
27529 


DteLTa 


0145 
e551 
e176 
el728 
0 2453 
23218 
23991 
©4756 
5490 
06175 
06795 
- 7333 
evttl 
21987 
«1705 
07173 


SAMPLES OF COMPUTER OUTPUT, 


DELTA-STAR 


0919 
20068 
©0135 
©0214 


NELTA-STAR 


©9018 
~0064 
0128 
©9202 
0282 
00344 
00442 
©0585 
0646 
-0701 
00747 
0779 
©0769 
©0706 
054? 


oo 





PROGRAM II 


THETA 








THETA 


-0014 
-0053 
-0107 
-0173 
0248 
- 0328 
oU4t1 
-0495 
-0578 
0656 
-0728 
-0£46 
- 0891 


.0927. 


9961 


2.161E 
1. 103€ 
2. 684E 





3-976E 


2-914E 
1.487E 
3-619E 


665 94E 


1.030E 
1Le-461E 
WAG 1) = 
2-438E 
22950t 
3.453E 
3.929E 
4. 360E 
4.731E 
5.-O028E 
5.240E 
Sie SELLE 


REYNOLDS NO. 


03 
04 


04 


REYNOLDS NO. 


03 
04 


Y-NORM/ 
0.22825E-01 


(0-53822E-01 


0.78683E-O1 
0.96202E-01 
0.10806E-00 
0-11424E-00 
0.11641E-00 
0.11685E-00_ 


oe eee ee § 


On ll7 2oE-00 
Os ll7T6Ze=00 
0-11794E-00 


(0.11821E-00 


0.11497E-00 
0-96544E-O1 


~ 0-69213E-01 


0.34408E-01 


Y-NORM/ 
0.-22821E-O01 
0.53803E-O1l 
0. /8642E-O1 


0.96138E-O1 


0.10797E-00 
0.11413E€-00 
0.11628E-00 
0.11670E-00 
0.11709E-00 
0.11744E-00 
0.11775E-00 
0.11L801E-00 
0.11477E-00 
0.96 362E-01 
0.69068E-O1 
0. 34343E-O1 





BOCY NUMBER= 
B8O0CY DEPTF=?2 
VELOCITY(KTS 


X( INCH) 
0.17 
0.87 
2-12 
rs} 
6-02 
8.53 
11.31 
14.25 
Ll7e 24 
20.19 
22-96 
25248 
27.65 
29.35 
30.63 
31.33 


BODY NUMGbER= 
BODY DEPTF=2 
VELOCITY(KTS 


X( INCH) 
0.19 
Os 97 
2-36 
4.30 
6.72 
9.53 

12.62 

15.90 

19.24 

22rd 

25462 

28.43 

30.85 

32. 79 

34.17 

34.96 


BODY NUMBER= 
BODY CEPTK=3 
VELOCITY(KTS 


X( INCH) 
0.22 
1.13 
2.75 
5.00 
7.61 

11.08 

14.69 

18.50 

22.39 

26. 20 

29.81 

33.08 

35-89 

38.14 

39.75 

40.67 


1 


)=0.75 


RADUTS 
O37 
0.88 
1.30 
1.60 
1.81 
le 96 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
1.63 
0.78 


3 


)=0.75 


KADUES 
0.40 
0.94 
leo 
1.67 
MST 
1.97 
2-00 
2-00 
2-60 
2-90 
2.09 
2.010 
1.94 
1.62 
1.15 
0.57 


fe) 


)=0.75 


RADUIS 
0.42 
1.00 
1.42 
1.72 
1-91 
2.00 
2.00 
2.00 
2.00 
2.00 
1.98 
1.85 
1.57 
1.18 
Ue %3 
3.28 


CORR. 
0.371 
0.880 
1.309 
1.621 
1.839 
1.998 
220%9 
2.054 
2-061 
2-067 
2.073 
2.078 
2-082 
2.C85 
1.697 
0.805 


CCRR. 
0.401 


CORR. 
0.422 
1.008 
1-436 
1.746 
1.946 
2-046 
2.055 
2-064 
2-072 
2.079 
2.065 
1.938 
1.657 
1.260 
0.795 
2316 


RADIUS 


RADIUS 


RADIUS 


DELTA 
0152 
01134 
~-1815 
©2576 
3379 
©4195 
~4997 
ooT6% 
2 648% 
7131 
e769 
8173 
08543 
8628 
07922 


DELTA 
0167 
«0610 
21234 
- 1980 
2808 
04559 
©5426 
26256 
«7030 
27730 
8340 
8826 
9051 


(© 8923 > 


8028 


DELTA 


20639 
21392 
°2228 
23155 
4130 
25113 
~ 6076 
e 70C0 
-/857 
° 8626 
29239 
9614 
©9655 
09154 
«7428 


CELTA=STAR 


~0067 
-0134 
-O2 1) 
©0295 
-0389 
20462 
20539 
~0610 
©0674 
-0729 
SORT 
-0816 
-0846 
-0818 
-0637 


DELTA-STAR 


-0020 
29073 
00146 
~0230 
29320 
0410 
0497 
-0579 
©0654 
0721 
0780 
~ 0830 
~0863 
20848 
2O7TT2 
~0586 


DELTA=STAR 


0023 
0082 
00163 
~O257 
©0356 
0456 
©0550 
0638 
0718 
«0790 
-0851 
0887 
~0882 
0824 
©0692 
©0443 


100 


THE Ta 
0015 
©0055 
©0113 
0182 
©0261 
0346 
©0433 
sOo72 
~0609 
~0692 
0767 
0835 
0892 
0937 
-0972 
© 1905 


THETA 
«0016 
« 0060 
0123 
0199 
20285 
«0378 
20474 
-0570 


«0665 


~ 0756 
-0839 
- 0912 
20975 
~ 1027 
«1069 
©1110 


THETA 


- 0018 
-0068 
-0139 
0225 
©0322 
«0426 
20535 
- 0645 
20752 


20854 


0948 
- 1033 
- 1106 
©1225 
21297 


6-626E 





2.932E 


4.585E 


6.497E 
Be 614E 
1.085F 
L231 3E 


mits > 386 


1.749E 
1.941E 
2-106E 
2.238E 
2e3356 
2.386E 


1.447E 
7.388E 
1.797E 
3.275E 
5-118E 
72258E 
9.612E 
1.-211E 


“1. 465E. 


1.715€ 
1.951E 
2e165E 
2-350E 
2-497E 
Ze 603E 
2. 663E 


1.676E 
8.606E 
2-094E 
3. 808E 
5. 948E 
8.439E 
Le1l19E 
1.409E 
1.705E 
1.995E 
2.270E 
2-519E 
22.734E 
2.905E 
3-027E 
3.-098E 


REYNOLDS NO. 
be 299e 03" 


03 


1.615E 04 


=| 


04 


REYNOLDS NO. 


REYNOLDS NO. 


03 


Y-NORM/ 
0.-23524E-01 
0.56149E-01 
0.82974E-01 


6. 16272E-00 


0.11656E-00 
0.12661E-00 
0.12967E-00 
0.13016E-00 


0. 13061E-00 





0.13101E-00 
0.13137E-00 
0.1L3167E-00 
O.13191E-00 
0.13210E-00 
0.10754E-00 


0.50993E-OL. 


Y-NORM/ 
0.22831E-O01 
0.53850E-01 
0.78741E-01 
0.96293E-O01l 
0.-10819E-00 
0O.11439E-00 
0. 11659E-00 
0.11 706E-00 
0. LI748E-00 
Ge. 1 7ATE-0e 
O.11820E-00 
0O.11849E-00 
0.-11524E-00 
0.69414E-01 © 


Y-NORM/ 
0. 20604E-01 
0.49250E-01 
0. 70191E-01 


“0. 85316E-01 


0.95092E-O1 
0.99979E-0l 
0.10044E-00 
0. LO087E-00 
0.10126E-0C0 
0.10161E-00 
0. 10093E-00 
0.94742F-01 
0.81L000E-01 
0.61602E-01 
0. 38848E-01 
0.-15430E-01 





BONY NUMBER=1 


CEPTH=2 
SPEED=0.75 


X(1) 
Om 73940E O00 
~0-94458F 00 
-0-86563E 00 
-0.-75540E OO 
=UeCl?Gs& OO 
-0.-45802E-60 
-~0-28160E-00 
=Ueo99013E-Ol 
Oo o013E-01 
0-28160E-00 
0-458C2E-00 
Oe61788E OO 
0O.75540E 00 
0-86563E O00 
0-94458E 00 
0.98940E 00 


BODY NUMBER=3 


CEPTH=2 
SPEED=0.75 


X(1) 
-0.98940E 00 
-0.94458F O00 
-0.86563E£ 00 
-0.75540E 00 
-0-61788E 00 
-0-45802E-CO 
-0.28160E-00 
-0.95013E-01 

0.95013E-01 
0-28160E-00 
0.45802E-00 
0-617828E 00 
0.75540E GO 
0.86563E 00 
0.94458E 00 
0.98940E 00 


TABLE G—-II — SAMPLES OF COMPUTER OUTPUT, PROGRAM III 


Y(i) 
0.23524E-Cl 
0-56149E-Cl 
0.829/74E-Cl 
Cel@7f2b-Co 
O-1LL656E-CO 
0.-126¢1E-CO 
Os IZ9E7E-CU 
0.13016F-CO 
O-13001E-CO 
O0.-13101E-CO 
O<1s157E-CO 
O<1310TE-CO 
Ono Le =CO0 
0.-13210E-CO 
O-10754E-CO 
GO. 309 73E-C) 


Y(1) 
0.27831E-Cl 
0.9385CE-Cl 
C.75741E-Cl 
0.964293E-Cl 
0O.1C819E-CO 
C.1L439E-CO 
C-11659E-CC 
O.L1706E-CO 
0O.11748E-CO 
Ce.lL17s87E-CO 
Calli zce—Co 
Ce-11849E-CO 
0.-11524E-CO 
C.96797TE-C]l 
0.69414E-Cl 
C.3449HE-Cl 


UCT) 
0.66402E 
0-98608E 
0. 74393E 
0. LO4SG63E 
Oe LOSSES 
O-10413E 
0.10222E 
0. 1LOLO6E 
0. LOOYCE 
O0-10099E 
O-1O136E 
0. 1LO200E 
O-10552ZE 


—O-12431E 


O-1L11424E 
0.55872E 


UCT) 
O.-701L15E 
0.100427F 
0.10542E 
0O.10473E 
O-1LO4LOE 
O-1C281E 
QO. LO1L66F 
0O-1L10088E 
O-LOOSILE 
O.10123E 
0.10062E 
0O.10435E 
O-LLO6SE 
O-1L0737€ 
0.-10298E 
0. 77644E 


GO 
OU 
OO 
Ol 
Ol 
Ol 
Ol 
Ol 
Ol 
01 
Ol 
Ol 
Ol 
Ol 
Ol 
00 


00 
Ol 
Ol 
Ol 
Ol 
Ol 
Ol 
01 
Ol 
Cl 
Ol 
Ol 
Ol 
Ol 
Ol 
CO 
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P( 1) 
0-55907E OO 
0.27646F-01 
0-43910E-00 
—-0.94643E-01 
=U. (231 1E-O1 
—Ce04595E—0 1 
-0.44932E-01 
=0. 2b251E=01 
=0.17995E-01 
-0.19858E-O1 
-0.27446E-01 
-0.40472E-01 
-0.11343E-00 
—0-.54533E 00 
-0. 30503E-00 

0-68783E OO 


P(1) 

0.-50839E OO 
=0. J9CIGCE-OZ 
=Ooll }25E-00 
=02-76953E-01 
=O 83725E—-01 
=O. 5.690 7E-C I 
-0. 33406E-01 
=O. L7OB9E-0] 
=O-16251F—-O1 
-O.24782E-01 
-0. 1644%8E-01 
-0.88884E-01 
=0.2259 176-00 
=O. 152792E-00 
-C.60430E-Ol 
0. 39714E-00 


PCC) 


O2927139E 
O.51912E 
O.525592E 
0.51722E 
0.51 756E 
0-51 738E 
O.5S1L7T99E 
0.51836E 
O.51341E 
0.-51838E 
0-51 826E 
0-51 806E 
0.51692E 
0-51020E 
0.51394E 
0452939 


PCC) 


0.-52660E 
0.51856E 
0.51696E 
O.51718E 
O.51739E 
0.51 780E 
O-SI8ILTE 
O.S51841E 
0.51844E 
0.51830E 
0-51643E 
O.57731E 
O-51518E 
O-51631E 
O-51775E 
0.52487E 


VEL(I)/ 
0.-49802E-00 
0.73956E 00 
0-56170E OO 
0-78469E O00 
O.-77682E OO 
O-78101E O00 
0.76666E O00 
OsTa17oe- Ue 
0.75672E OO 
0.75741E 00 
0.76022E 00 
0.76503E OO 
O0.79140E O00 
0.93233E 00 
0.-85678E 00 
0.41904E-00 


VEL(I)/ 
0-52586E O00 
O./75318E O00 
0.-79062E O00 
O-78551E OO 
O.-78077E OO 
O.77104E 00 
0.-76242E O00 
0. T5060E 700 
O0.75606E O00 
O-75924E 00 
0. 75614E O00 
0.-/78262E O00 
0-83016E 00 
0.805376 00 
0-77233E 00 
0.58233E O00 


BODY 
DEP TH=1 
SPEED=1.12 


X{( I) 
-0.98940E CO 
-0.-94458E 00 
-0.86563E GO 
-0.75540E 00 
Sumeol?8s8E OO 
-0-45802E-00 
-0.28160E-00 
mur J5013E-01 

Oe 959013E-Ol 
0.28160E-00 
0.45802E-00 
O.61788E O00 
0.75540E 00 
WeeeoG3E 00 
0-94458E 00 
0.-98940E 00 


BODY NUMBERK=3 


GEPTH=2 
SPEED=1.51 


x(1) 
-0.98940E 00 
-0.94458E 00 
-0.86563E 00 
-0.75540E 00 
-0.61788E 00 
~0.45802E-00 
~0.28160E-00 
-0.95013E-Ol 

0.95013E-01 
0.28160E-00 
0 .45802E-00 
0.61788E 00 
0.75540E 00 
0.86563E 00 
0-94458E 00 
0.98940E 00 


BODY NUMBEK=5 


DEPTH=3 


SPEED=0.75 


X{( I) 
-~0.98940E CO 
—-0.94458E 00 
—-0.-860563E 00 
-0.75540E GO 
~0.-61788E 00 
-0.45802E-00 
=O. 20 1606-01) 
=O, ID 01L3E=0 | 

0.95013E-01 
Oe 25 LO0E=00 
Oa 7b 04e—-00 
0O.61788E OO 
0-75540E 00 
0.86563E O00 
912.94458E 00 
0.98940E 00 


NUMBEK = 3 


Y(t) 
O0-Z22502Z5E-U1 
Oe o3ecce-O0l 
Us7GO55E—-Cl 
O226202E-Cl 
0-L0806E-CO 
O.-L1L424E-CO 
0.-11641E-CO 
Oe PlGe5E-00 
O0.-11725E-00 
O.117&62E-CO 
Os bli 74e-CC 
O-LIdZ2&lE-CO 
O-LL497E-CC 
0.96544E-01 
0.69215E-Cl 
0.34408E-Cl 


Vi) 
Oe22G21C=Cl 
0.53803E-O1l 
0.78642E-01 
GO.26136F-O01 
OeLC 19 7E-CO 
0.11413E-00 
0.-1L1628E-00 
OslT6G710E-00 
Oc L1L709eE-CO 
O.11744E-CO 
OS Lit 7 5E=CO 
0O.-1LIL8G1LE-CO 
Os 1LES77E-CO 
0.96362E-C1l 
0.-69068E-01 
0.34343E-01 


Y(1) 
MeacUOUGE-O1 
0.49250E-O01l 
Oe 7 Ol TE-O 
Oso S316E-—Cl 
0. J509Z7E-C I 
Ovo oT IE—C 1 
0.1L10044E-CO 
0. LCO87TE-—CO 
eLOlLZoe-CO 
GeLtOlLéle—co 
0. 2009 2E-CO 
0.94742E-Cl 
0.-S8LOO0OE-Cl 
USoloOoZe—- Cl 
0.38848EF-Ol 
Del s430e—Cl 


UCT) 
O27 O1STE 
0.10043E 
0.10542E 
0. 10473E 
O0.-LO409E 
O-L0280E 
O-1LCLO5E 
0. LOO8TE 
0O-10079C 
0. LOOY6E 
O-1LOL4S9E 
O-1L0399E 
O-LILO64E 
O-LO732E 
0.-1L0289E 
0.776056 


Ul) 
0. 70143E 
0.10043E 
0.10542E 
0O-1L0473E 
0O-LO409E 
0.1023d0E 
O-lLOL65E 
0. LOO8GE 
O-LOO79E 
O. TO095E 
O-L0148E 
0-.10398F 
O.-11L062E 
Celie Jt 
O-LU284E 
O. 77 580E 


Ut.) 
Os 703 71E 
0.102856 
0.-10469t 
0.104355 
O21 039 E 
O~-1LO267E 
0.-10093E 
O-LOO7LE 
0. LOO76E 
0.-10133E 
OC. LOS09E 
0. 10432E 
0.-10346E 
O-1LOL35E 
OS I2G1CE 
O.771 38E 


P(1) 

O.50816E O00 
=O. 861726E-02 
=O LEN ZOE-CO 
= (9.969) 1LE =O.) 
=e GS46TE—O01 
eG. 2OSO9E-OL 
=Ueu 2156-01 
=O. b7421E-01 
—Uslo789E—O] 
—OnlI25I9E—O! 
aie 27 OL Ou 
-0.-41449E-01 
=0.22402E-00 
=O. Fo169E=—00 
=O. 3075 1E-0) 

Qed 7 15E- 00 


P(1) 

0.50799E 00 
-0.87155E-02 
-0.11127E-00 
-0.96905E-01 
-0.83397E-O1 
~0.56824E-01 
-0.33229E-01 
-0.17329E-O1 
-0.15799E-01 
-0.19173E-01 
-0.29892E-01 
-0.81180E-O1 
-0.22358E-00 
-0.15105E-00 
-0.57650E-O1 

0. 39813E-00 


P(1) 
0.50479E OO 
-0.57780E-0O1 
-0.96035E-O01 
-0.88918E-0O1 
-0.71388E-01 
-0.54101E-Ol 
-0.18611E-O1 
-0.14236E-0O1l 
-0.15212€-01 
-0.26690E-01 
-C.62655E-01 
-0. 88368F-Ol 
-0. 70363E-O01 
-0.27248E-Ol 
0.82038E-0l 
0.40498E-00 
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PCC(I) 


O69 L35E 
0-93340€E 
0.92984E 
0<,93033E 
0.93080E 
0293 L72E 
0.93254E 
O27 933098 
O.93315E 
C.93303E 
0.93266E 
GeI3ZOBKBTE 
0-92592E 
0.92843E 
0.93 166E 
0.94750E 


BCE 
0.55072E 
0.51814E 
O-5LL6E7E 
0.51258E 
0-51343E 
0O-51510E 
0.-51659E 
0.51 760E 
Oss T69E 
0.51 74HE 
O.51680E 
0.51357E 
0.50459E 
0.50916E 
0.51505E 
0.54379E 


PCC(I1) 
0.-19462E 
O. LO SaTE 
OwLSSZ27E 
0.-18538E 
O. 18566E 
OO. LS 5935 
0. 18648E 
0.-18655E 
0. 18653E 
0.186355 
Osi So TIE 
0.18539E 
0O.-1L8567EF 
0.18634E 
0.-18804E 
O- 19307E 


— 


VEL Cli 
0. /78547E OO 
O-11248E Ol 
O-11807E OL 
O-LL?730E OL 
O-1L1658E OL 
O-11514€ Ol 
O.-11385E Ol 
O.11297E Ol 
O-11288E Ol 
O-LL3Z07E Ol 
O-1LL367E Ol 
O-11L647E Ol 
O.-12391E Ol 
O.12019E Ol 
O-LLS24E Ol 
0.-86917E OO 


VEEL 
Os 1OS59ZE 20 
O-15166E Ol 
Os Goa eR. ] a en | 
O-L5815E Ol 
Ota tiie Ot 
0.15523E Ol 
0.15349E Ol 
0.152306 01 
Oe oZ19E Ol 
O-15244E Ol 
0.15324E Ol 
0O-15701t Ol 
O-16703E Ol 
0.16200E Ol 
Os LOSZIE O01 
O.L1L/15€ O2 


VEL(I)/ 
0.5277, 00 
O- T71SGE 00 
Ge 7b519%e 00 
0.78263E OO 
C.7To3SLE OO 
0.77 002E -00 
0.75095t OO 
O../Soaze 00 
O01 To 560E- 00 
0.75394E 90 
O.-77314E OO 
0.-78244E 00 
O-77594E OOD 
0.76015E O00 
O. T1L8>8e) 00 
0.57453E CO 
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